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Abstract 

We study the critical behavior (CB) of all period p-tuplings (p = 2, 3, 4, . . .) 
in (A^ = 2,3,4, . . .) symmetrically coupled one-dimensional maps. We first 
investigate the CB for the N = 2 case of two coupled maps, using a renor- 
malization method. Three (five) kinds of fixed points of the renormalization 
transformation and their relevant "coupling eigenvalues" associated with cou- 
pling perturbations are found in the case of even (odd) p. We next study the 
CB for the linear- and nonlinear-coupling cases (a coupling is called linear or 
nonlinear according to its leading term), and confirm the renormalization re- 
sults. Both the structure of the critical set (set of the critical points) and the 
CB vary according as the coupling is linear or nonlinear. Finally, the results 
of the two coupled maps are extended to many coupled maps with > 3, in 

which the CB depends on the range of coupling. 
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I. INTRODUCTION 



Universal scaling behavior of period p-tuplings {p = 2, 3, 4, . . .) has been found in a one- 
parameter family fA{x) of one-dimensional (ID) unimodal maps with a quadratic maximum. 
As the nonlinearity parameter A increases, a stable fixed point undergoes the cascade of 
period-doubling bifurcations accumulating at a finite parameter value Aoq. The period- 
doubling sequence corresponding to the MSS (Metropohs, Stein, and Stein |]T|) sequence 
i?*" (for details of the MSS sequences and the (*)-composition rule, see Refs. exhibits 
an asymptotic scaling behavior 0. 

What happens beyond the period-doubling accumulation point Aoo is interesting from the 
viewpoint of chaos. The parameter interval between Aoo and the final boundary-crisis point 
Ac beyond which no periodic or chaotic attractors can be found within the unimodality 
interval is called the "chaotic" regime. Within this region, besides the period-doubling 
sequence, there are many other sequences of periodic orbits exhibiting their own scaling 
behavior. In particular, every primary pattern P (that cannot be decomposed using the 
(*)-operation) leads to an MSS sequence P*". For example, P = RL leads to a period- 
tripling sequence, P = RL^ to a period-quadrupling sequence, and the three different P = 
RLR^, RL^R, and RL^ to three different period 5-tupling sequences. Thus there exist 
infinitely many higher period p-tupling (p = 3,4, . . .) sequences inside the chaotic regime. 
Unlike the period-doubling sequence, stability regions of periodic orbits in the higher period 
]9-tupling sequences are not adjacent on the parameter axis, because they are born by their 
own tangent bifurcations. The asymptotic scaling behavior of these (disconnected) higher 
period p-tupling sequences characterized by the parameter and orbital scaling factors, 6 and 
a, vary depending on the primary pattern P P,^|-|TO[]. 

In this paper we consider (A^ = 2,3,4,...) symmetrically coupled ID maps, which 
may be used as models of coupled nonlinear oscillators such as Josephson-junction arrays 
or chemically reacting cells, and so on We are interested in the critical behavior (CB) 
of period p-tuplings (p = 2, 3, . . .) in the coupled ID maps. The period-doubling case with 
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p = 2 was previously studied |T2|JT^. Here we study the critical scaling behavior of all the 
other higher period p-tuplings (p = 3, 4, . . .). 

Using a renormalization method, we first investigate the critical behavior for the N = 2 
case of two coupled maps in Sec. 0. In the case of even (odd) p, we find three (five) kinds 
of fixed points of a renormalization transformation and their relevant coupling eigenvalues 
(CE's) associated with coupling perturbations. A short account of the renormalization 
result has been already published [Q. We next consider two kinds of couplings, linear- 
and nonlinear-coupling coupling is called linear or nonlinear according to its leading 

term. As examples of the linear- and nonlinear-coupling cases, we study the linearly- and 
dissipatively-coupled maps, respectively in Sec. |T|, and confirm the renormalization results. 
The structure of the critical set (set of the critical points) varies depending on the nature 
of coupling. In the linearly-coupled case, an infinite number of critical line segments and 
the zero-coupling critical point, at which the ID maps become uncoupled, constitute the 
critical set, while in the dissipatively-coupled case, the critical set consists of only one critical 
line segment, one end of which is the zero-coupling critical point. The CB also depends on 
the position on the critical set. For even (odd) p, three (four) kinds of fixed points govern 
the CB for the linearly-coupled case, whereas only two (three) fixed points govern the CB 
for the dissipatively-coupled case. 

In Sec. 1^ we extend the results of the two coupled maps to many coupled maps with 
A^ > 3. It is found that the critical scaling behavior depends on the range of coupling. In 
the extreme long-range case of global coupling, in which each ID map is coupled to all the 
other ID maps with equal strength, both the structure of the critical set and the CB are 
the same as those for the two-coupled case, independently of A^. However, for the cases of 
nonglobal couplings of shorter range, a significant change in the structure of the critical set 
may or may not occur according as the coupling is linear or not. For the case of a linear 
nonglobal coupling, only the zero-coupling critical point is left in the parameter plane, which 
is in contrast to the global-coupling case. On the other hand, for the case of a nonlinear 
nonglobal coupling, one critical line segment still remains, as in the globally-coupled case. 
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Finally, a summary is given in Sec. 0. 



II. RENORMALIZATION ANALYSIS OF TWO COUPLED MAPS 

In this section we first discuss stability of periodic orbits in two coupled ID maps, and 
then study the CB associated with period p-tuplings {p = 2,3,4,.. .) using a renormalization 
method. Three (five) kinds of fixed points of a renormalization operator and their relevant 
CE's with moduli larger than unity are found in the case of even (odd) p. 



A. Stability of periodic orbits in two coupled maps 

We consider a map T consisting of two symmetrically coupled ID maps, 

_ I xt+i = F{xt, yt) = f{xt) + g{xt, yt), 
1 : < (2.1) 

[ yt+i = F{yt,xt) = f{yt) + g{yt,xt), 

where t denotes a discrete time, f{x) is a ID unimodal map with a quadratic maximum at 
X = 0, and g{x, y) is a coupling function. The uncoupled ID map / satisfies a normalization 
condition 

/(O) = 1, (2.2) 

and the coupling function g obeys a condition 

g{x,x) = for any re. (2.3) 

The two-coupled map ( p.l|) is invariant under the exchange of coordinates such that 
X ^ y. The set of all points which are invariant under the exchange of coordinates forms a 
symmetry line y = x. An orbit is called a(n) (in-phase) synchronous orbit if it lies on the 
symmetry line, i.e., it satisfies 

Xt = yt for all t. (2.4) 
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Otherwise, it is called an (out-of-phase) asynchronous orbit. Here we study only the syn- 
chronous orbits, which can be easily found from the uncoupled ID map, Xt+i = f{xt), 
because of the condition (|2.3[). 

Stability analysis of a periodic orbit in the two-coupled map T can be conveniently 
carried out in a set of new coordinates (X, Y), defined by 



_ (x + y) y _ {x- y) 



(2.5) 



Here X and Y correspond to the synchronous and asynchronous modes of the orbit, respec- 
tively. For example, for a synchronous orbit X = x and F = 0, whereas for an asynchronous 
orbit r 7^ 0. 

In order to study the stability of a synchronous orbit with period g, we consider an 
infinitesimal perturbation (5X, 5Y) to the orbit. Here 6X and 5Y correspond to the 
synchronous-mode and asynchronous-mode perturbations, respectively. Linearizing the q- 
times iterated map T'^ (expressed in terms of the new coordinates) at an orbit point, we 
obtain a linearized map. 



5Y, 



t+q 



J 



(2.6) 



where the Jacobian matrix J (= DT'^) of T"? is given by the g-product of the Jacobian matrix 
DT of T along the orbit: 



9-1 

J=Yl DT{xt 



n 

t=0 







V 







(2.7) 



f'ixt) - 2Gixt) 

Here the prime denotes a derivative with respect to x, and G{x) = dg{x,y)/dy \y=x', here- 
after, G{x) will be referred to as the "reduced coupling function" of g{x,y). Note that 6X 
and 6Y become decoupled for the case of a synchronous orbit (i.e., J has a diagonalized 
form) . 
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The eigenvalues of J, called the stability multipliers of the orbit, are then given by: 



The two stability multipliers, Aq and Ai, determine the stability of the synchronous orbit 
against the synchronous- mode and asynchronos-mode perturbations, respectively. Hereafter, 
they will be called the synchronous and asynchronous stability multipliers, respectively. 
Note also that the synchronous stability multiplier Aq is just the stability multiplier of the 
uncoupled ID map, and the coupling affects only the asynchronous stability multiplier Ai. 

A synchronous orbit is stable when it is stable against both the synchronous-mode and 
asynchronous-mode perturbations, i.e., the moduli of both stability multipliers are less than 
unity (|Aj| < 1 for i = 0,1). Hence, the stable region of a synchronous orbit in the parameter 
plane is bounded by the synchronous and asynchronous bifurcation lines determined by the 
equations Aj = ±1 for i = 0,1, as will be seen in Sec. fT|. When the Aq = 1 (—1) line 
is crossed, the synchronous orbit loses its stability via synchronous saddle-node (period- 
doubling) bifurcation. On the other hand, when the Ai = 1 (—1) line is crossed, it becomes 
unstable via asynchronous pitchfork (period-doubling) bifurcation. Some brief explanations 
on the bifurcations are given below. 

In the case of a synchronous saddle-node bifurcation, the synchronous orbit collides with 
an unstable synchronous orbit with the same period, and then they disappear, like the tan- 
gent bifurcation in the ID maps. On the other hand, there are two types of supercritical and 
subcritical bifurcations for each case of the pitchfork and period-doubling bifurcations. In 
the supercritical case of the synchronous (asynchronous) pitchfork and period-doubling bi- 
furcations, the synchronous orbit loses its stability, and gives rise to the birth of a pair of new 
stable synchronous (asynchronous) orbits with the same period and a new stable synchronous 
(asynchronous) period-doubled orbit, respectively. However, in the subcritical case of the 
synchronous (asynchronous) pitchfork and period-doubling bifurcations, the synchronous 
orbit becomes unstable by absorbing a pair of unstable synchronous (asynchronous) orbits 
with the same period and an unstable synchronous (asynchronous) period-doubled orbit. 




(2.8) 



t=o 



t=0 
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respectively. (For more details on bifurcations on 2D dissipative maps, refer to ||15|| .) 



B. Renormalization analysis 

We now consider the period p-tupling renormalization transformation A/", which is com- 
posed of the p-times iterating (T*^p^) and rescaling (B) operators: 

Af{T) = BT^P^B-\ (2.9) 

Here the rescaling operator B is: 



B 



la 0^ 

y a J 



(2.10) 



because we consider only synchronous orbits. 

Applying the renormalization operator M to the coupled map ( |2.1| ) n times, we obtain 
the n-times renormalized map T„ of the form, 

(2.11) 

Vt+i = Fn{yt,xt) = fniyt) + gnivuxt). 
Here fn and Qn are the uncoupled and coupling parts of the n-times renormalized function 
Fn, respectively. They satisfy the following recurrence equations: 

/„+i(x)= afi^\-), (2.12) 
a 

9n+^{x,y) = aF^\^, 1) - afl^\-), (2.13) 
a a a 

where 

fi'\^) = fn{fi'-'\x)), (2.14) 

F(^)(x,i/)= Fr,{Flf-'\x,y),Flf-'\y,x)), (2.15) 

and the rescaling factor is chosen to preserve the normalization condition /„+i(0) = 1, i.e.. 



The recurrence relations ( p.l2|) and (|2.13|) define a renormalization operator TZ of transform- 
ing a pair of functions (/, g): 



(f \ 

Jn+l 
9n+l I 



n 



( f \ 

Jn 



(2.17) 



A map Tc with the nonhnearity and couphng parameters set to their critical values is 
called a critical map: 



(2.18) 



xt+i = Fc{xt,yt) = fc{xt) +gc{xt,yt), 
Vt+i = F^iVuXt) = fcivt) + 9ciyt,xt). 

A critical map is attracted to a fixed map T* under iterations of the renormalization trans- 
formation J\f: 



xt+1 = F*{xt,yt) = f*ixt) + g*{xt,yt), 
yt+i = F*{yt,xt) = f*{yt) + g*{yt,xt). 

Here {f*,g*) is a fixed point of the renormalization operator TZ with a = 1/ f*^P~^\l): 



(2.19) 









= n 









(2.20) 



This fixed-point equation can be solved row by row consecutively. Note that f*{x) is just the 
fixed function in the ID map case, which varies depending on p Only the equation 

for the coupling fixed function g*{x,y) is therefore left to be solved. One trivial solution is 
g*{x,y) = 0. In this zero-coupling case, the fixed map ( p.l9|) consists of two uncoupled ID 
fixed maps, which is associated with the CB at the zero-coupling critical point. 

However, it is not easy to directly find coupling fixed functions other than the zero- 
coupling fixed function g*{x, y) = 0. We therefore introduce a tractable recurrence equation 
for a reduced coupling function G{x) = dg{x,y)/dy\y=x. Differentiating the recurrence 
equation ( p.l3| ) for g{x,y) with respect to y and setting y = x, we obtain a recurrence 
equation for G{x) [jl6|: 
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a 



a a a 

+f'r''\-)Gn{fr'\-)). 



(2.21) 



where F^^{x) = dFjf\x,y) /dy\y=x- Then Eqs. ( p.l2| ) and ( |2.21| ) define a "reduced renor- 
mahzation operator" TZ of transforming a pair of functions (/, G) : 

/ 

V 



fn+1 


\ 

= n 


( f 

J n 


\ 

(2.22) 


Gn+1 


1 


\Gn 


/ 


of n, 


which satisfy 


(A 


/ 




(2.23) 


= n 


\^ J 


\ 


G^] 





Here /* is just the ID fixed function and G* is the reduced couphng fixed function of g*, 
i.e., G*{x) = dg*{x,y)/dy \y=x. 

We first consider the lowest two cases with p = 2,3. For p = 2, Eq. (|2.21| ) becomes 



Gn+i{,x) 



f(2) 



.X, 



a 



Gni-Mifni-)) - 2GMn{-))] 



a 



a 



a 



+ /n(-)G.(/„(-)). 

a a 



(2.24) 



We note that ^^^(f ) in Eq. (g^) becomes 



p(2) 



,X , 



a 



for Gn{x) = 0, 



n,2 \ 



a 



for GM = ^fni^) 



7(2), 
n,2 \ 



a 



.X, 



a 



(-)-l] 
z a 



a 



for G„(x) = -[/:(x)-l], 

^[/f (-)-!] 
/ a 

for G„(x) = ^[/;(x) + l], 
for Gn{x) = f'n{x). 
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(2.25a) 
(2.25b) 

(2.25c) 

(2.25d) 
(2.25e) 



Differentiating the ID fixed-function equation for f*{x) with respect to x in the period 
p-tupling case, we have 



p-i 



/*'(x) = r(p)'(-) = nr'(r^^H-))- 



(2.26) 



Using Eqs. ( p.25|) and (|2.26|) , we obtain three solutions for the reduced couphng fixed function 



G*(x) 



G*{x) = 0, 
G*ix) = ^r'ix), 
G*ix) = \[r'ix)-1]. 



(2.27a) 
(2.27b) 
(2.27c) 



Thus we find three fixed points (/*, G*) of TZ for the period-doubling case. 

We next consider the period-tripling case with p = 3. The recurrence equation for G 
becomes 



a„+i(x) = Fi5(-) 

a 



= i^i2(-)[/;(/f(-))-2G„(/J^)(:^))] 
a a a 

a a 

We also note that ^^^(f ) in Eq. (gl^ becomes [see Eq. for 
^(3). a;. 



n,2\ 



forG„('x) = 0, 



a 



a 



7(3), 
n,2 1 



^ X ^ 



^/f'(-) for G.(x) = ^/;:(x), 
Z a z 

(-)-!] 



1 



f,S(-) 

q; 



i^S(-)=/f'(-) forG„(x) = /:(x). 



for G.(x) = -[/;(a;)-l], 

^[/f (^) + l] 

for G„(x) = ^[/:(x) + l], 

, X , 



Using Eqs. (|2.26| ) and (|2.29|) , we get five solutions for G*{x): 



(2.28) 



(2.29a) 
(2.29b) 

(2.29c) 

(2.29d) 
(2.29e) 



10 



G*{x 
G*{x 
G*{x 
G*{x 
G*(x 



0, 

^[r'(^)-i], 



1 



[r'(^) + i], 



2 



(2.30a) 
(2.30b) 
(2.30c) 
(2.30d) 
(2.30e) 



Thus we find five kinds of fixed points (/*, G*) for the period-triphng case. 

We now consider the general period p-tupling case. By induction, we obtain the solutions 
for G*{x) as follows. Assume that F,lf2{-) for the case of period p-tupling becomes 



7(f) ( 
n,2\ 



.X , 



a 



a 



a 



for GJx) = 0, 

for GJx) 



1 r{p)' /^-^ 

2^" 



2 a 
for GJx) = ■ 



;[/::(^)-i], 



(2.31a) 
(2.31b) 

(2.31c) 



for even p 
+ for odd p 



a 



for Gn{x 




(2.31d) 



for even p 

f(py{^) for odd p 
for GM = r^{x). 



(2.31e) 



For the case of the next period (p + l)-tupling, the recurrence equation (|2.21|) for G becomes 



a 

= i^S(-)[/;(/i^n-))-2G„(/„(^)(-))] 
a a a 



(2.32) 



It follows from the assumption ( |2.31| ) that i^i^2^^'*(f ) for the case of period (p + l)-tupling 
also satisfies Eq. (|2.31|) . Note that the lowest even and odd cases, p = 2 and 3, satisfy the 
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assumption (|2.31|) [see Eqs. ( |2.25|) and ( p.29| )]. Consequently, the equation ( p.31j ) is valid for 



all p. Using Eqs. ( |2.26| ) and ( p.31D , we obtain three (five) solutions for G*{x) in the case of 
even (odd) p\ 

G*(x)=0 for all p, (2.33a) 

G*{x) = ^r'{x) for all p, (2.33b) 

G*{x) = i[/*'(x) - 1] for all p, (2.33c) 

G*{x) = ^[f*\x) + 1] for odd p, (2.33d) 

G*{x) = f*\x) for odd p. (2.33e) 

Thus we find three (five) kinds of fixed points (/*, G*) of TZ for the case of even (odd) p. 
In case of a critical map ( p.l8|) , the synchronous and asynchronous stability multipliers 



Ao,n and Ai_„ of the synchronous orbits are given by [see Eq. (pl8|)] 

Ao,n = ii fci^t), \i,n = lilf'cixt) " 2G,{xt)], (2.34) 

i=0 t=0 

where Gc{x) is the reduced coupling function of gdx), i.e., Gc{x) = dgc/dy\y=x- As n ^ oo, 
they converge to their limit values Ag and A*, called the critical synchronous and asyn- 
chronous stability multipliers, respectively: 

Ag = lim Ao,„, A* = lim Ai,„. (2.35) 



n— >oo n— >oo 



Here Aq is just the critical stability multiplier A* for the ID case, and the coupling affects 
only \\. The values of A* depend on the fixed points, as shown below. 

The invariance of a fixed map T* under iterations of the period p-tupling renormalization 
transformation M implies that, if T* has a periodic point (x, y) with period p", then B^^{x, y) 
is a periodic point of T* with period p^'^^ . Since rescaling leaves the stability multipliers 
unaffected, all synchronous orbits of period {n = 0,1,2,...) have the same stability 
multipliers, which are just the critical stability multipliers. Consequently, they have the 
values of the stability multipliers of the fixed point of the fixed map T*: 
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(2.36) 



where x is the fixed point of the ID fixed function [i.e., x = f*{x)], and Aq is just the critical 
stabihty multipher A* of the uncoupled ID map. Note that XI depends on the reduced 
coupling fixed function G*{x): 



XI = X* for G*{x) = 0, 

XI = foTG*{x) = ^r'{x), 

XI = 1 forG*(x) = ^[r'(a:)-l], 

XI = -1 forG*(x) = ^[r'(x) + l], 

A* = -A* for G*{x) = f*'{x), 



(2.37a) 
(2.37b) 
(2.37c) 
(2.37d) 
(2.37e) 



where the cases (|2.37a|) - (|2.37c| ) exist for all p, but the cases (|2.37d|) - (|2.37e|) exist only for 
odd p. 

Consider a pair of functions {f*,G). Here G{x) is not necessarily a reduced coupling 
fixed function G*{x). When fn{x) = f*{x) and Gn{x) = G{x), the function F.l:^l{^) of 
Eq. i ^m^j will be denoted by Fi^^(f ). We now examine the evolution of a pair of functions 
(/* + /i, G + $) close to (/*, G) under the reduced renormalization operator TZ. Linearizing 
TZ at {f*,G), we obtain a linearized operator C of transforming a pair of infinitesimal 
perturbations (/i, $): 



where 



''-n.+l 



n+1 



c 



\*»/ 





C2 



I 



(2.38) 



hn+l{x) = [Clhn]ix) 

a 



a 



y ^ H ) ] linear 

a 



a a a 



^n+l{x) = [C2^n]ix) + [Cshn]ix) 



(2.39) 
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n,2 ' 



a 



[£2$n](x) =[/*'(/ 



a 



a 



2 V ) ] linear ; 

a 



a a 



a 



)), 



(2.40) 



(2.41) 



a 



'a' 



a 



a 



)) 



a 



a 



(p-i)' 



q; 



a 



'a' 



(2.42) 



Here the variations Sfjf^^) and 5-Fi^2(^) ^.re introduced as the hnear terms (denoted by 
_ /*(p)(|)]H,ear and [^i^^^f ) - F^^^ ( f )] linear ) m h and $ of the deviations of /(^')(f ) 
and Fi5(f ) from f*^P\^) and ^^^^(f ), respectively. 

When G{x) is a reduced fixed couphng function G*{x) of Eq. (|2.33|) , the operator C of 
Eq. (p.38|) becomes a hnearized transformation of 7^ at a fixed point {f*,G*). A pair of 
perturbations {h*, $*) is then called an eigenperturbation with eigenvalue u, if it satisfies 



(2.43) 



The reducibility of C into a semiblock form implies that to determine the eigenvalues of C 
it is sufficient to work independently in each of h{x) subspace and $(x) subspace. That is, 
one can find eigenvalues of £i and £2 separately and then they give the whole spectrum of 
C. 

We first solve the eigenvalue equation for £1, i.e., 









' h* 


\ 






= c 

















uh*ix) = [Cih*]{x). 



(2.44) 



Note that this is just the eigenvalue equation for the ID map case. It has been shown 
that there exists only one relevant eigenvalue 5, associated with scaling of the nonlinearity 
parameter, whose values vary depending on p P,^|-[TO[|. However, note that although the 
eigenvalue 6 of £1 is also an eigenvalue of £, {h, 0) is not an eigenperturbation of C unless 
£3 is a null operator. 
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We next consider a perturbation of the form (0, $) having only the couphng part. In 
this case (0, $*) can be an eigenperturbation of £, only if satisfies 



(2.45) 



Eigenvalues associated with coupling perturbations are called CE's. 

For the case of coupling perturbation (0, $) to (/*, G), C2 of Eq. ( p.41|) becomes: 



a 



a 



a' 



i-*(p-i)/^ 

a' 



))• 



(2.46) 



We first consider the lowest two cases with p = 2, 3. For the period-doubling case, 6F2^\^) 
becomes 



6Fi'\^ 
a 



SF 



(2), 



a 

(2),X 



sFr'i 



a 

SF.PC- 

a' 



r'(r(-))$(-) + r'(^)$(/*(^)) for g{x) = o, 

a a a a 

forG(x) = ^r'(x), 

H-) + Hn-)) ioTG{x) = hr'{x)-i], 

a a I 

-$(-) - $(/*(-)) for G(x) = \\r\x) + 1], 
a a 2 

-/*'(/*(-))$(-) - /*'(-)$(/*(-)) for G{x) = r\x) 
a a a a 



(2.47a) 
(2.47b) 
(2.47c) 
(2.47d) 
(2.47e) 



Note that only G"s for the first three cases (|2.47a|) - (|2.47c| ) are the reduced coupling fixed 
functions G*'s. Using Eqs. and (|^, we also find 5F!f\^) for p = 3, 



(3) 



5Fr 



5fP 



6fP 



SF 



(3) 



-) = y /*W'(^)$(/*W(^))/*(2-0'(^*(^+l)(^)) _ Q 



a 



i=Q 



-) = forG(x) = V(x), 
a 2 

-) = E-^-ir^^H-)) for G{x) = lir'ix) - 1], 
-) = E<^>(r^^n-)) for G{x) = lir'ix) + 1], 



^) = Er^^^'(^)'^(r^^n^))r^'"^^'(r^^^'H^)) for g(x) = r'(x) 



i=0 



(2.48a) 
(2.48b) 
(2.48c) 

(2.48d) 

(2.48e) 
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where f*^^\x) = x. For this period-triphng case, all G's in Eqs. ( p.48a| ) - ( p.48e|) are the 
reduced coupling fixed functions G*'s. 

We obtain SF2^\^) for the general period p-tupling case by induction. Suppose that 
5^2^"* (^) for the period p-tupling case is 



5Ff (-) 

a 



p~i 



for G{x) = 0, 



i=0 



a 



a 



a 



6R 



ip), 



a 



a 



for = -r (a;), 

forG(x) = l[r'(x)-l], 



p-i 



for even p 



a 



a 



X, 



i=0 

p-1 



.X 



y $(/*«(-)) for odd p 
p-i 



for G(x) = -[r'(x) + l], 



(2.49a) 
(2.49b) 
(2.49c) 

(2.49d) 



i=0 

p-i 



^ /*W'(-)$(/*«(-))/*(f-i-^)'(/*(^+i)(-)) for odd p 



for G(x) = f*\x). 



a 



a 



(2.49e) 



For the case of the next period (p + l)-tupling, 5F^^^^ of Eq. (|2.46|) becomes 



5Fi^+i) = \r\f<P\-))~2GU*^P\-))]5F!f\-) 
a a a 



+ [/HP)'(n)_2Fi^)(-)]$(/*w(-)). 

a a a 

It follows from Eqs. (p.31|) and (|2.49|) that 5F!f'^^'^ for the period {p + l)-tupling case also 
satisfies Eq. ( |2.49|) . Note that the lowest even and odd cases, p = 2 and 3, satisfy the 
assumption (|2.49|) [see Eqs. ( |2.47|) and ( p. 4^ )]. The equation ( |2.49|) is therefore valid for all 
p. 

We now find CE's associated with coupling perturbations. For the zero-coupling case of 



(2.50) 



G*{x) = 0, the CE equation ( |2.45|) becomes 
z/$*(x) = [C2^*]{x) = 6Fi 



ip) I ' 



a 



a a a 



(2.51) 
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Using the fact that /*'(0) = 0, it can be easily shown that when x = 0, the equation (p.51|) 
becomes 

i^^*{0) = ['[ir'{f<'\0))]^*{0). (2.52) 

i=l 



Letting x — in Eq. ( p^.26D , we also have 



nVnr^^Ho)) = lim|^ = a. (2.53) 

1=1 J \a) 

Then Eq. (|2.52| ) reduces to 

zy$*(0) = a$*(0). (2.54) 

There are two cases. For the case $*(0) ^ 0, we have the first CE 

Pi = a. (2.55) 

The eigenfunction $J(a;) with CE ui is of the form, 

$*(x) = 1 + flix + azx^ H , (2.56) 

where a^'s {i = 1,2, . . .) are some constants. For the other case $*(0) = 0, it is found that 
f*'{x) is an eigenfunction for the CE equation (|2.51|) . When ^*(x) = f*'(x), the equation 
( p.51| ) becomes 

i^r\x) = pf<^y{-) = pr\x). (2.57) 

a 

We therefore have the second relevant CE 

1^2=?: (2.58) 

with reduced coupling eigenfunction 

m^) = r'i^)- (2.59) 

Note that $2(^) no constant term, while $*(x) has a constant term. Thus we find two 
relevant CE's, 1/1=0 and 1^2 = p, for the zero-coupling case. 
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The nth image $„ of a general reduced couphng perturbation $ under the hnear trans- 
formation C2 has the form, 

~ aivl^l{x) + a2V2^*2{x) for large n, (2.60) 

because the irrelevant part of $„ becomes negligibly small for large n. Here ai and 02 are 
relevant components. 

A coupling is called linear or nonlinear according to its leading term. In the case of a 
linear coupling, in which the coupling perturbation (p{x,y) has a leading linear term, its 
reduced coupling function $(x) has a leading constant term, and hence $(0) 7^ 0. However, 
for any other case of nonlinear coupling with a leading nonlinear term, its reduced coupling 
function has no constant term, and hence $(0) = 0. Note that the relevant component ai 
becomes zero for the nonlinear-coupling case, while it is nonzero for the linear-coupling case. 
Consequently, the CB associated with coupling perturbations is governed by two relevant 
CE's z/i = a and 1^2 = P for the linear-coupling case, but by only one CE 1^2 = P for the 
nonlinear-coupling case, which will be confirmed in Sec. [TTH 

We now consider the cases of two (four) other reduced coupling fixed functions (j*(x)'s 



for even (odd) p [see Eq. ( p.33|) ]. They are associated with CB at critical points other 
than the zero-coupling critical point, as will be seen in Sec. |T|. Since SF2^\-) = for 



G*{x) = |/*'(x), C2 becomes a null operator. Hence there exist no CE's, and the CB is 
essentially the same as that for the ID case. When G*{x) = |[/*'(x) — 1], the CE equation 
( p.45| ) becomes 



z/$*(x) = \C2^*]{x) = SF^''\- 

a' 

p-i 



^ $*(/*«(-)). (2.61) 



i=0 " 



There exists a relevant CE 
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(2.62) 



when $*(a;) is a nonzero constant function, i.e., = 6 (6 is a nonzero constant). For 

odd p there are two additional reduced couphng fixed functions, G*{x) = |[/*'(a;) + 1] and 
G*{x) = f*'{x). The CE equation (|2^) for G*{x) = l[f*'{x) + 1] is just that of Eq. ( p^H ). 
Therefore it has the same CE u = p as that for the case G*{x) = l[f*'{x)-l]. However, the 



critical asynchronous stability multipliers XI for the two cases are different [see eq. p.37| )]. 
When G*{x) = f*'{x), the CE equation (|2^ ) is the same as that for the case G*{x) = 0. 
Hence it also has two relevant CE's ui = a and 1^2 = P- However, \l for this case is different 
from that for the case G*{x) = 0, as can be seen in Eq. (|2.371 ). The results of relevant CE's, 
along with those of the critical asynchronous stability multipliers, are listed in Table |. 



III. LINEAR AND NONLINEAR COUPLINGS 

We choose f{x) = 1 — Ax"^ as the uncoupled ID map and consider two kinds of cou- 
plings, linear and nonlinear couplings. As examples of the linear- and nonlinear-coupling 
cases, we study the linearly- and dissipatively-coupled maps, respectively, and confirm the 
renormalization results. The structure of the critical set (set of the critical points) varies 
depending on the kind of coupling. In the linearly-coupled case, an infinite number of criti- 
cal line segments, together with the zero-coupling critical point, constitute the critical set, 
whereas in the dissipative case, the critical set consists of the only one critical line segment, 
one end of which is the zero-coupling critical point. The CB also depends on the position 
on the critical set. For even (odd) p, three (four) kinds of fixed points govern the CB for 
the linearly-coupled case, while only two (three) kinds of fixed points govern the CB for the 
dissipative case. 



A. Linearly coupled maps 

We numerically study the CB of period p-tuplings in two linearly-coupled ID maps with 
the coupling function 
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9ix,y) = ^(y - x)^ (3-1) 

where c is the couphng parameter. The critical scahng behavior depends on whether p is 
even or odd. 

As an example of odd period p-tuplings, we take the period-tripling case {p = 3) and 
study its CB. The stability diagrams of synchronous orbits with period g = 3" (n = 0, 1, 2, 3) 
are shown in Figs. Q-]^. As noted in Sec. |11 A| , the stable region of each synchronous orbit 
of level n (period 3") in the parameter plane is bounded by the four synchronous and 
asynchronous bifurcation curves determined by the equations Aj^„ = ±1 for i = 0,1. Since 
each synchronous orbit of level n is born by its own synchronous saddle-node bifurcation 
(which occurs for Ao,n = 1), a sequence of stability regions with increasing n is not connected, 
unlike the period-doubling case [0. We now examine the treelike structure of stability 
regions. Figure |l| shows the stability regions of synchronous orbits of the lowest two levels 
(i.e., n = and 1). The synchronous orbit with period g = 1 is stable in some quadrilateral- 
shape region containing the c = line segment. However, the stability region of the next 
level n = 1 consists of three quadrilateral-shape areas. They can be regarded as "daughter" 
quadrilaterals of the "mother" quadrilateral of level 0. That is, it may be thought that they 
branch off from the mother quadrilateral. 

We next consider the stability regions of higher levels in Fig. ^ The branchings occur 
from the central one containing the c = line segment and its nearest neighboring one (i.e., 
the right one) among the three quadrilaterals of level 1 [see Figs, ^(a) and|^(b)], while there 
is no branching from the left one [see Fig. ^(c)]. This rule governs the treelike structure of 
the stability regions. That is, for each level n branchings occur only from two quadrilateral- 
shape areas, the central one containing the c = line segment and its nearest neighboring 
one. However, an infinite number of successive branchings occur only for the case of the 
central quadrilateral including the c = line segment [see Fig. 0(a)]. For the case of the 
nearest neighboring quadrilateral, branching occurs only once, and after that, successive 
quadrilaterals of higher levels pile up without any further branchings [e.g., see Fig. @(b)]. 
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For the cases of quadrilaterals other than the central and its nearest neighboring ones, 
successive quadrilaterals of higher levels pile up without any branchings [e.g., see Fig. 0(c)]. 

A sequence of stability regions with increasing period is called a "period-tripling route" , 
like the period-doubling case [|13|. There are two kinds of period-tripling route. The sequence 
of the quadrilateral-shape areas containing the c = line segment converges to the zero- 
couphng point c = on the A = A^^) line, where A^^ (= 1.786 440 255 563 639 354 534 447 . . .) 
is the accumulation point of the period-tripling sequence for the ID case. It will be re- 
ferred to as the route. On the other hand, a sequence of quadrilaterals which piles 
up without successive branchings converges to a critical line. For example, the left- 
most one is the line joining two points q (= —3.590 291636 032 974400 442...) and 
(= -3.482 633 674 606 564177473...) on the A = A^^^ line [see Fig. |(c)]. This kind of 
route will be called an L3 route. Note that there are infinitely many L3 routes, while the 
Z3 route converging to the zero-coupling critical point {Aoo, 0) is unique. Hence an infinite 
number of critical line segments, together with the zero-coupling critical point, constitute 
the critical set. 

We now study the critical scaling behavior on the critical set. First, consider the case of 
the Z3 route ending at the zero-coupling critical point. The CB for this zero-coupling case 
is governed by the zero-coupling fixed map (|2.19| ) with g*{x,y) = 0, which has two relevant 
CE's ui = a {= -9.277341 . . .) and 1^2 = 3 (see Table |). 

The quadrilateral-shape areas of the levels n = 1, 2, 3 in the Z3 route are shown in Fig. |^. 
Note that Figs. |5(a), ^(b), and ^(c) nearly coincide except for small numerical differences. 



It is therefore expected that the height and width hn and Wn fl^ of the quadrilateral-shape 
area of level n may geometrically contract, 

hn ~ <5~", Wn ~ a"" for large n, (3.2) 

where S = 55.247026 .... Thus the scaling factors of the nonlinearity and coupling param- 
eters become 6 and a, respectively, which will be shown explicitly below. 

We follow the synchronous orbits of period g = 3" up to level n = 10 in the route, 
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and obtain a self-similar sequence of parameters {An, c„), at which each orbit of level n has 
some given stability multipliers (Aq, Ai) (e.g., Aq = —1 and Ai = 1). Then the sequence 
{{An,Cn)} converges geometrically to the zero-coupling critical point {A^,0). In order to 
see the convergence of each of the two scalar sequences {An} and jcnj, define 5„ = 
and Un = , where AAn = An — An-\ and Ac^ = c„ — c„_i. Then they converge to their 
limit values 6 and /i, respectively, as shown in Table |T|. Hence the two sequences {^n} and 
{c„} obey one-term scaling laws asymptotically: 

AA„ ~ 5"", Ac„ ~ /i"" for large n, (3.3) 

where b = 55.247... and fi = —9.277... . Note that the nonlinearity-parameter scaling 
factor 6 is just that for the ID case, and the value of the coupling-parameter scaling factor 
fi is close to that of the first CE z/i (= a). 

In order to take into account the effect of the second relevant CE i>2 (= 3) on the scaling 
of the sequence {Ac„}, we extend the simple one-term scaling law ( |3.3[ ) to a two-term scaling 
law ITHl: 



Ac„ ~ Ci/i^ " + C2IJ.2 " for large n, (3.4) 

where \^2\ > \f^i\, and Ci and C2 are some constants. This is a kind of multiple scaling law 
0. Eq. O gives 

ACn = SiAc„+i - S2Ac„+2, (3.5) 

where Si = fii + /i2 and S2 = /ii/i2- Then /xi and /i2 are solutions of the following quadratic 
equation, 

/x^ - Sin + S2 = 0. (3.6) 

To evaluate /xi and /X2, we first obtain si and S2 from Ac„'s using Eq. (^.51): 
_ Ac„Ac„+i - Ac„_iAc„+2 _ Ac^ - Ac„+iAc„_i 

Sl — 5 , S2 — -7 K 7 7 . I'J- ' j 

Ac;+i - Ac„Ac„+2 Ac;+i - Ac„Ac„+2 
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Note that Eqs. (|3.4|) -(|3r7| ) hold only for large n. In fact the values of Sj's and /ij's (i = 1, 2) 
depend on the level n. Therefore we explicitly denote Sj's and /ij's by Sj^^'s and iii^n^i 
respectively. Then each of them converges to a constant as n — > oo: 

lim Si^n = Si, lim /ij^^ = /ij, i = 1, 2. (3.8) 



Three sequences {/ii,n}, {Ai2,n}, and {jJ^^nl (^"i,-"} are shown in Table |T|. The second 
column shows rapid convergence of to its limit value /ii (= —9.277341 . . .), which is 
close to the renormalization result of the first relevant CE vi (= a). (Its convergence to a is 
faster than that for the case of the above one-term scaling law.) From the third and fourth 
columns, we also find that the second scaling factor /i2 is given by a product of two relevant 
CE's vi and 1^2, 

/i2 = -, (3.9) 

where vi = a and z/2 = 3. It has been known that every scaling factor in the multiple- 
scaling expansion of a parameter is expressed by a product of the eigenvalues of a linearized 
renormalization operator . 



We also study the effect of CE's on the asynchronous stability multipliers of synchronous 
periodic orbits. Consider the two-coupled map (|2.1| ) with /(x) = fc{x) and g{x,y) = 
eip{x, y). Here fc{x) is the ID critical map with the nonlinearity parameter set to its critical 
value A = and e is an infinitesimal coupling parameter. The map for e = is just the 
critical map at the zero-coupling critical point consisting of two uncoupled ID critical 
maps fc- It is attracted to the zero-coupling fixed map ( |2.19| ) with F*{x,y) = f*{x) under 
iterations of the renormalization transformation A/" of Eq. ( p.9| ). Hence the reduced coupling 
function G{x) [= = ed{p{x,y)/dy\y=x] corresponds to an infinitesimal reduced coupling 
perturbation to the reduced coupling fixed function G*{x) = 0. 

In the period p-tupling case, the stability multipliers Ao,n and Ai,„ of the ^"-periodic orbit 
are the same as those of the fixed point of the n times renormalized map A/'"'(T), which are 
given by 
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Ao,n = fLi^n), Ai,„ = f^{xn) - 2G„(x„) ^ f^{xn) - 2e$(£„). (3.10) 

Here (/„, be the nth image of (/c, G) under the reduced renormahzation transformation 
TZ, Xn is just the fixed point of fn{x) [i.e., x„ = fn{£n)] and converges to the fixed point 
X of the ID fixed map f*{x) as n —* oo. The first stabihty multipher Ao,n converges to 
the ID critical stabihty muhipher A* = /*'(£) as n — >■ oo. For the period-triphng case, 
A* = -1.872 705 929. . . . Since - [QG]{x) = e$„(x) [<l>„(x) is given in Eq. (g^)], 

the asynchronous stabihty muhipher has the form 

~ A* + e [eiz/^" + eaz/^] for large n, (3.11) 

where ei = —2ai^l{x) and 62 = —2a2f*'{x). Therefore the slope Sn of Ai^„ at the zero- 
coupling critical point (e = 0) is 

Q _ ^Ai,, 
oe 



^eii/" + e2Z/2 for large n. (3.12) 

£=0 

Here the coefficients ei and 62 depend on the initial reduced function $(a;), because a„'s 
are determined only by $(x). Note that the coefficient ei is zero for the nonlinear-coupling 
case, whereas it is nonzero for the linear-coupling case. Hence the growth of 5"^ for large n 
is governed by the two relevant CE's ui = a and 1^2 = P for the linear-coupling case, but 
only by the second relevant CE 1^2 = p for the nonlinear-coupling case. 

Figure ^ shows three plots of Ai^„(y4^\ c) versus c for n = 2, 3, and 4. The slope Sn of 
Xi^n at the zero-coupling critical point increases with n, and obeys well the two-term scaling 
law, 

Sn ~ (ii^i + d2r2 for large n, (3.13) 
where di and ^2 are some constants and |ri| > \r2\. This equation gives 

Sn+2 = tlSn+1 — t2Sn, (3-14) 

where ti = ti + r2 and t2 = rir2- As in the scaling for the coupling parameter, we first 
obtain ti and t2 of level n from S'„'s: 
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, _ Sn+lSn ~ Sn+2Sn-l , _ *S'„S'„_|_2 , s 

Then the scahng factors ri „ and r2^n of level n are given by the roots of the quadratic 
equation, — ti,„r„ + t2,n = 0. They are listed in Table and converge to constants ri 
(= ui) and r2 (= 1^2) as n ^ 00, whose accuracies are higher than those of the coupling- 
parameter scaling factors. 

We next consider the cases of L3 routes, each of which converges to a critical line segment. 
In each L3 route, there are two kinds of self-similar sequences of parameters {A^, c„), at which 
each orbit of level n has some given stability multipliers; the one convergs to the left end 
point of the critical line segment and the other convergs to the right end point. As an 
example, consider the leftmost L3 route [see Fig. H(c)], which converges to the critical line 
segment with two ends (A^^ , ci) and {A^^ ^ Cr). We follow, in the leftmost L3 route, two self- 
similar sequences of parameters, one converging to the left end and the other converging to 
the right end. In both cases, the sequence {An} converges geometrically to its accumulation 
value A'^^ as in the case of the route, 

AA„~(5~" for large n, (3.16) 

where Ay4„ = An — An-i and 5 = 55.247. . . . The sequences {cn} for both cases also obey 
the one-term scaling law, 

Ac„ ~ for large n, (3.17) 

where Ac„ = c„ — c„,_i. The convergence of the scaling factor /i„ of level n to its limit value 
/i (= 3) is shown in Table [V|. Note that the value of fi is different from that (/i = a) at the 
zero-coupling critical point. Although the scaling factors of the coupling parameter at both 
ends are the same, the critical asynchronous stability multipliers A*'s of Eq. ( p. 351) at both 



ends have different values. The convergence of the sequence {Ai,n} to its limit value A* is 
also shown in Table ^ At the left (right) end, A^ = 1 (— 1). Comparing the values of /i and 
A^ with those of the CE u and A* listed in Table |, we find that the CB near the left end is 
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governed by the fixed point (/*, G*) of TZ witli G*{x) = — 1], wliereas tliat near tlie 

riglit end is governded by tlie fixed point with G*{x) = ^[f*'{x) + 1]. 

Figure |] shows the behavior of the asynchronous stabihty multipher Ai^„(A^\ c) near the 
leftmost critical line segment. The slopes S^s of Ai „ at both ends obey well the one-term 
scaling law, 

Sn ~ z^"" for large n, (3.18) 

where z/ = 3. For any fixed value of c inside the critical line segment, Ai^„ converges to 
zero as n ^ oo. That is, all the interior points are critical points with A* = 0. Hence the 
CB inside the critical line segment becomes the same as that of the ID map, which will be 
discussed in more details below. This kind of ID-like CB is governed by the fixed point with 
G*{x) = which has no relevant CE's (see Table |). 

For the case of a synchronous orbit in two linearly-coupled ID maps, its two Lyapunov 
exponents are given by 

1 m—l 

aoiA)= hm -J2\n\f'ixt)\, (3.19) 
aM,c) = hm - ^ \n\f'ixt)-c\. (3.20) 

m^oo rn — 

Here ao(ai) is the synchronous (asynchronous) Lyapunov exponents characterizing the mean 
exponential rate of divergence of nearby orbits along (across) the symmetry line y = x. 
Note that the synchronous Lyapunov exponent ctq is just that of the uncoupled ID map 
and the coupling affects only the asynchronous Lyapunov exponent ai. In order to see 
the phase dynamics near the critical line segment in more details, we fix the value of the 
nonlinearity parameter A = A'^^ and obtain the asynchronous Lyapunov exponent ai of 
the synchronous orbit by varying the coupling parameter c, which is shown in Fig. ^. (Note 
that the synchronous Lyapunov exponent ctq of the synchronous quasiperiodic orbit becomes 
zero on the A = A^^ line.) Inside the critical line segment (q < c < Cr), the synchronous 
quasiperiodic orbit on the y = x symmetry line becomes a synchronous attractor with ai < 0. 
Since the dynamics on the synchronous attractor is the same as that for the uncoupled ID 
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case, the critical maps at interior points exhibit essentially ID-like CB. However, as the 
coupling parameter c passes through c; or c^,, the asynchronous Lyapunov exponent cti of 
the synchronous quasiperiodic orbit increases from zero, and hence the coupling leads to 
desynchronization of the interacting systems. Thus the synchronous quasiperiodic orbit 
ceases to be an attractor outside the critical line segment, and new asynchronous attractors 
appear. This is illustrated in Fig. |^. 

We also study the critical scaling behavior of the asynchronous Lyapunov exponent ai 
near both ends for the case A = . As shown in Fig. |, the asynchronous Lyapunov 
exponent cxi varies linearly with respect to c near both ends, i.e., (Xi ~ £, e = c — c* 
(c* = Ci or Cr). The critical scaling behavior of ai near both ends are obtained from the 
same CE u = 3 oi the fixed points with G*{x) = |[1 ± Consider a map with nonzero 

e (but with A = A^^^) near both ends. It is then transformed into a new one of the same 
form, but with a renormalized parameter e' under a renormalization transformation. Here 
the parameter e obeys a scaling law, 

£' = ue = 3e. (3.21) 

Then the asynchronous Lyapunov exponent (Xi satisfies the homogeneity relation, 

aiie')=3a,{e). (3.22) 

This leads to the scaling relation, 

al(£)~£^ (3.23) 

with exponent 

= ln3/lnz/ = 1. (3.24) 

We now study the CB of period quadruplings (p = 4), as an example of even period-p 
tuplings. The stability diagrams of synchronous orbits of period g = 4" (ra = 0, 1,2,3) are 
shown in Figs. |9|-pi?|. The treelike structure of stability regions is similar to that for the period- 
tripling case. As shown in Fig. four "daughter" quadrilaterals of level 1 branch off from 
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its "mother" quadrilateral of level 0. However, an infinite number of successive branchings 



occur only for the case of the central quadrilateral including the c = line [see Fig. |TO|(a)]. 
For the cases of quadrilaterals other than the central one, successive quadrilaterals of higher 
levels pile up without any branchings [see Figs, ^(h) and 0(c)]. 

Like the period-tripling case, there are two kinds of period-quadrupling route which is 
a sequence of stability regions with increasing period. The sequence of the quadrilaterals 
containing the c = line segment, called the route, converges to the zero-coupling critical 
point {A^^\ 0), where A^^) (= 1.942 704 354 755 467 972 167 178 . . .) is the accumulation point 
of the period-quadrupling sequence for the ID case. On the other hand, a sequence of quadri- 
laterals which piles up without branchings, called an L4 route, converges to a critical line. For 
example, the leftmost one is the line joining two points q) and (A^^ , c^) [see Fig. p!0|(c)], 
where q = -3.888 058 931 772 634 488 .. . and = -3.877 063 178 096 222 051 ... . Note also 
that there are infinitely many L4 routes. Hence an infinite number of critical line segments, 
together with the zero-coupling critical point, constitute the critical set. 

We first consider the case of the Z4 route ending at the zero-coupling critical point. The 
critical scaling behavior for this zero-coupling case is governed by the zero-coupling fixed 
map with two relevant CE's Ui = a {= —38.819 074 . . .) and 1^2 = 4. 

Figure |ll| shows the quadrilaterals of the levels n = 1, 2, 3 in the Z4 route. Note that 



Figs. |TT](a), |TT|(b), and [TT|(c) nearly coincide. The height and width /i„ and w„ of the 
quadrilateral of level n is therefore expected to geometrically contract, 

hn ~ 5""", Wn ~ for large n, (3.25) 

where 6 = 981.594 976 .... Thus the scaling factors of the nonlinearity and coupling param- 
eters become 6 and a, respectively, which will be shown in more details below. 

We follow, in the Z4 route, the synchronous orbits of period g = 4" up to level n = 7, 
and obtain a self-similar sequence of parameters (y4„,Cn), at which each orbit of level n 
has some given stability multipliers (Aq, Ai). Then the sequence {{An, c^)} converges to the 
zero-coupling critical point (A^\0). Like the period-tripling case, the two sequences {An} 
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and {c„} obey well one-term scaling laws, 

Av4„ ~ (5"", Ac„ ~ /i"" for large n, (3.26) 

where AAn = y4„ — An-i and Ac„ = c„ — c„_i. The scaling factors 5„ and /i„ of level n 
converge to their limit values 5 (= 981.594 . . .) and ^ (= —38.819 . . .), respectively, as shown 
in Table VI. Note that 5 is just the nonlinearity-parameter scaling factor for the ID case, 
and the value of the coupling-parameter scaling factor /i agrees well with that of the first CE 
vi (= a). In fact, the second CE V2 (= 4) also has an effect on the scaling of the sequence 
{Ac„}. Thus it obeys well the two-term scaling law, 

Ac„ ~ Ci/i-" + Cs/i^" for large n, (3.27) 

where |/i2| > and Ci and C2 are some constants. As shown in Table VII, the two scaling 
factors /ii and ^2 are given by 

/^i = ^1, /i2 = — • (3.28) 

We also study the effect of CE's on the asynchronous stability multipliers of synchronous 
orbits. Figure |T2] shows three plots of Ai_„(yl^\ c) versus c for n = 1, 2, 3. The slope Sn of 
Xi^n at the zero-coupling critical point obeys well the two-term scaling law, 

Sn ~ dir'l + dara for large n, (3.29) 



where di and ^2 are some constants and |ri| > \r2\- As shown in Table |VIII| , the scaling 



factors Ti n and r2,„ of level n converge to constants ri (= Vi) and r2 (= ^2)- 

We next consider the cases of L4 routes, each of which ends at a critical line segment. 



As an example, consider the leftmost L4 route [see Fig. |TO|(c)], in which we follow two self- 
similar sequences of parameters, one converging to the left end (A^^q) of the critical line 
segment and the other one converging to the right end (A^^c^). As in the case of the Z4 
route, the sequence converges geometrically to its limit value with the ID scaling 

factor 5. The sequence {c„} also obeys the one-term scaling law, 
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Ac„ ~ /i "for large n, (3.30) 

where Ac^ = c„ — c„--i. The convergence of the scahng factor /i„ of level n to its limit value 
is shown in Table Note that the scaling factors /i at both ends are the same, i.e. 



= 4. Moreover, the critical asynchronous stability multipliers A]^'s at both ends are also 
the same, i.e. A* = 1, as shown in Table This is in contrast to the period-tripling case 
where at both ends are different. We also compare the values of yU and \\ with those 
of the CE V and \\ listed in Table H and find that the CB at both ends is governed by the 
same fixed point (/*,(?*) of Iz with G*{x) = — 1], unlike the period-tripling case. 

Figure |13| shows the behavior of the asynchronous stability multiplier Ai.„(A^\c) near 
the leftmost critical line segment. The growth of the slope S'„ of Ai „ at both ends is governed 
by the CE v (=4), i.e., 

Sn ~ z/" for large n. (3.31) 

However, for any fixed value of c inside the critical line egment, Ai „ converges to zero as 
n — i> oo. Thus, all the interior points become critical ones with A* = 0. Consequently, the 
CB inside the critical line segment becomes the same as that of the ID map, as will be seen 
below. This kind of ID-like CB is governed by the fixed point with G*{x) = ^/*'(x), which 
has no relevant CE's (see Table |). 

In order to see the phase dynamics near the critical line segment in more details, we fix 
the value of the nonlinearity parameter A = and obtain the asynchronous Lyapunov 
exponent ai [see Eq. ( |3.20|) j of the synchronous orbit by varying the coupling parameter c. 



which is shown in Fig. [T^. Inside the critical line segment (q < c < Cj.), the synchronous 
quasiperiodic orbit on the y = x symmetry line becomes a synchronous attractor with 
o"i < 0, as shown in Fig. [l^(a). Note that the dynamics on the synchronous attractor is 
the same as that for the uncoupled ID case. Hence the critical maps inside the critical line 
segment exhibit ID-like CB. However, as the coupling parameter c passes through q or Cr, 
the asynchronous Lyapunov exponent ai of the synchronous quasiperiodic orbit increases 
from zero, and hence the coupling leads to desynchronization of the interacting systems. 

30 



Thus the synchronous quasiperiodic orbit ceases to be an attractor outside the critical hne 
segment, and new asynchronous attractors appear, as shown in Figs. |l^(b) and p!5|(c). 

We also study the scaling behavior of the asynchronous Lyapunov exponent ai near both 
ends of the critical line segment. As shown in Fig. |T6|, the asynchronous Lyapunov exponent 
varies linearly with respect to c near both ends, i.e., ai ~ e, e = c — c* (c* = q or c^). As 
in the period-tripling case, using the scaling theory, one can also obtain the scaling relation 
of (Ti, 0"! ~ e'' with exponent 77 = ln4/ In z/ = 1, where u (= 4) is the CE of the fixed point 
with G*{x) = lir'ix) - 1]. 

Finally, we briefly summarize the results for the linear-coupling case. The critical set 
consists of the zero-coupling critical point and an infinite number of critical line segments. 
The CB at the zero-coupling critical point is governed by the zero-coupling fixed map with 
two CE's Ui (= a) and z/2 (= p) for all period p-tupling cases. However, the CB near both 
ends of each critical line segment depends on whether p is even or odd. In the case of odd p, 
the CB at one end is governed by a fixed point (/*, G*) of IZ with G*{x) = |[1 — f*'ix)] and 
that at the other end by another fixed point with G*{x) = |[1 + f*'{x)]. These two fixed 
points have only one CE v = p. On the other hand, in the case of even p, the CB at both 
ends is governed by the same fixed point with G*{x) = \[f*'{x) — 1]. Inside the critical line 
segment, the CB is the same as that of the ID map for all period p-tupling cases. This kind 
of ID-like CB is governed by the fixed point with G*{x) = ^f*'{x), which has no relevant 
CE's. Consequently, for even (odd) p, three (four) kinds of fixed points govern the CB for 
the linearly-coupled case. 

B. Dissipatively coupled maps 

As an example of the nonlinear-coupling case, we consider two dissipatively-coupled ID 
maps with the coupling function 

9ix,y) = ^[f{y)-f{x)], (3.32) 

and study the CB of the period triplings (p = 3) and period quadruplings (p = 4). 
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Figures 0(a) and |T7|(b) show the stabihty diagrams of synchronous orbits with period 
q = {n = 1,2,3) for p = 3 and 4, respectively. As previously shown, each stability 
region of level n (period p") in the parameter plane is bounded by four bifurcation curves 
determined by Aj^„ = ±1 for i = 0,1. An infinite sequence of such stability regions, called the 
"period p-tupling route" , converges to a critical line joining two ends ( A^-* , C;*) and {A^^^ , c*) , 
where A^^^ is the accumulation point of the period p-tupling sequence for the ID case, = 
and — 2, as shown below. Hence only one critical line segment constitutes the critical set 
for the dissipative case, unlike the linearly-coupled case. 

Consider two dissipatively-coupled ID maps on the line A = A^^} in the parameter plane, 
in which case the reduced coupling function of the couphng function (|3.32|) is given by 

G{x) = ^/^(x), (3.33) 

where fc is the ID critical map with the nonlinearity parameter set to its critical value. By 
successive applications of the renormalization operator TZ to {fc,G), we have 

Ux) = a/S(-), G„(x) = ^/;(x), (3.34) 
Cn = cj^_i - 34_i + 3c„_i for p = 3, (3.35) 
Cn = -Cn-i + 4c^_i - 6c^_i + 4c„_i for p = 4, (3.36) 

where fo{x) = fc{x), Go{x) = G{x), and cq = c. Here converges to the ID fixed map 

The fixed points of the recurrence equations (|3.35|) and ( p.36| ) for c are denoted by solid 
circles in Fig. For p = 3, there are three fixed points, 

c* = 0,1,2, (3.37) 

while only two fixed points, 

c* = 0, 1 , (3.38) 

exist for p = A. Stability of a fixed point c* is determined by its stability multiplier A given 
by A = dcn/dcn-i\c*- Note that the fixed point at c* = 1 is superstable (A = 0). The basin 
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of attraction to the superstable fixed point becomes tlie open interval (0,2), because any 
initial c inside the interval < c < 2 converges to c* = 1. For the period-tripling case, both 
ends cf = and c* = 2 of the interval are unstable fixed points with A = 3, and all points 
outside the interval diverge to the plus or minus infinity. However, for the quadrupling case, 
only the left end cj = 0, which is also the image of the right end c* = 2 under the recurrence 
relation (|3.36| ), is an unstable fixed point with A = 4, and all points outside the interval 



diverge to the minus infinity. Thus the line segment connecting two end points c^* = and 
c* = 2 becomes the critical line for the dissipative-coupling case. 

All critical maps inside the critical line segment are attracted to the fixed maps with the 
same reduced couphng function G*{x) = These fixed maps have no relevant CE's, 

because the fixed point c* = 1 is superstable. Hence the critical maps at interior points 
exhibit essentially ID-like CB. In the period-tripling case, the critical map at the left end 
cf is attracted to the zero-coupling fixed map, 

T* : Xi+i = r{xt), yt+i = rivt), (3.39) 

with the reduced coupling function G*{x) = 0, and the critical map at the right end c* to 
another fixed map, 

T* : Xi+i = rivt), yt+, = /*(xO, (3.40) 

with the reduced coupling fixed function G*{x) = f*' (x). However, in the period-quadrupling 
case, the critical maps at both ends are attracted to the zero-coupling fixed map (|3.39|) . Note 
that the two fixed maps of Eqs. (|3.39|) and (|3.40|) have the same relevant CE's ui = a and 
1^2 = P (see Table 1). However, for this dissipatively-coupled case, the CB near both ends are 
governed only by the second CE 1^2 = p (i.e., the first relevant component ai in Eq. (p.60|) 



becomes zero), which can be easily understood from the fact that the fixed points c* = 0, 2 
are unstable ones with stability multiplier X = p. 

We study the critical scaling behavior associated with coupling near the critical line 
segment and confirm the renormalization results. For the dissipative-coupling case, the 
stability multipliers of synchronous orbits with period g = p" {p = 3, 4) become 
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q-l 

Ao,„(A) = n f'i^t), c) = (1 - cyXo,n. (3.41) 

i=0 

Let the pair of stability multipliers of the synchronous orbit of level n (period p") at a point 
{A,c) be (Ao,Ai). Then there exists a "conjugate point" {A, —c + 2), at which the pair of 
stability multipliers becomes (Aq, — Ai) and (Aq, Ai) for p = 3 and 4, respectively. For c = 1, 
Ai = and the two conjugate points become degenerate. 

Like the linearly-coupled case, we follow the synchronous orbits of level n in the period 
p-tupling route, and obtain a self-similar sequence of parameters {An,Cn), at which each 
orbit of level n has some given stability multipliers (Aq, Ai). Without loss of generality, we 
choose Ao = — 1. Then one can find a pair of mutually conjugate sequences. One sequence 
{{An,Cn)} can be obtained by fixing — 1 < Ai < 0, which converges to the zero-coupling 
critical point {A^,0) as follows: 

A^ - A^P^ ~ r", c„ = 1 - (-Ai)P"" ~ - ln(-Ai)p-" for large n, (3.42) 

where A^^} and S are just those for the ID case. Note that the coupling-parameter scaling 
factor is just the second relevant CE 1/2 (= p) of the zero-coupling fixed map ( 3.39 ). The 



"conjugate sequence" can also be obtained by following ^"-periodic orbits with Ai_„ = — Ai 
and Ai for p = 3 and 4, respectively. This kind of conjugate sequence {(y4„,c„)} converges 
to the right-end critical point (A^\ 2) (i.e., the conjugate point of the zero-coupling critical 
point) as follows: 

An - ~ 5-'^, c„ - 2 = -1 + (-Ai)P"" ~ ln(-Ai)p-" for large n. (3.43) 

Note also that the asymptotic scaling of the coupling-parameter sequence is governed only 
by the second CE 1/2 = p of the "conjugate" fixed map ( pj.40| ) for p = 3 and of the zero- 
coupling fixed map ( |3.39| ) for p = 4. Hence the critical scaling behavior of the coupling 
parameter near both ends become the same. 

Figures |l9](a) and |l^(b) show three plots of Ai^„(y4^\c) versus c for = 1,2,3 in the 
period-tripling {p = 3) and period-quadrupling [p = 4) cases, respectively. The critical 
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asynchronous stability multipliers A*'s of Eq. (|2.35|) at both ends of the critical line segment 



can be easily obtained from Eq. ( p.41| ). For p = 3, = A*(— A*) at the left (right) end (A* 



is the critical stability multiplier for the ID case), while for p = 4, A^ = A* at both ends. 
The slopes S^s of Ai^„ at both ends also obey well the same one-term scaling law, 



Sn 



dc 



p"' for large n, (3.44) 



where c* = c* or c*. Hence the growth of Sn for large n is governed only by the second CE 
(=p)- 

All interior points of the critical line segment become critical ones with A* = 0, because 
at any interior point Ai^„ converges to zero as n ^ oo. Hence the CB inside the critical line 
segment becomes the same as that of the ID map, as will be seen below. This kind of ID-like 
CB is governed by the fixed maps with the same reduced coupling function G*' (x) = |/*'(x), 
which have no relevant CE's, as shown above. 

For the dissipatively-coupled case, the asynchronous Lyapunov exponent ai is given by 

ai(A,c) =ao(A) +ln|l-c|, (3.45) 

where ctq is the ID Lyapunov exponent of Eq. ( |3.19| ). In order to see the phase dynamics 
near the critical line segment, we fix the value of the nonlinearity parameter A = A';^. Then, 
the asynchronous Lyapunov exponent becomes ai{A^P^ , c) = In |1 — c|, because ao{A^) = 0. 
Inside the critical line segment (0 < c < 2), the synchronous quasiperiodic orbit on the y = x 
symmetry line becomes a synchronous attractor with ai < [see Figs. 0(a) and |l^(a)]. Note 
that the dynamics on the synchronous attractor is the same as that for the uncoupled ID 
case. Hence the critical maps inside the critical line segment exhibit ID-like CB. However, 
as the coupling parameter c passes through or c*, the asynchronous Lyapunov exponent 
0"! of the synchronous quasiperiodic orbit increases from zero, and hence the coupling leads 
to desynchronization of the interacting systems. Thus the synchronous quasiperiodic orbit 
ceases to be an attractor outside the critical line segment, and new asynchronous attractors 
appear, as shown in Fig. |2^. 
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We also study the scaling behavior of the asynchronous Lyapunov exponent cxi near 
both ends of the critical line segment. The asynchronous Lyapunov exponent near both 
ends becomes a\ ~ e, where e = c — c* (c* = or c*). As in the linearly-coupled case, using 
the scaling theory, one can also obtain the scaling relation of ai, o\ ~ with exponent 
r] = Inp/ In = 1, where z/2 (= p) is the second CE of the fixed maps ( |3.39| ) and ( ^.4(j| ). 



IV. EXTENSION TO MANY COUPLED MAPS 

In this section we study the CB of period p-tuplings in (A^ > 3) coupled ID maps, 
in which the coupling extends to the A'th [1 < K < ^i^Y^) even (odd) A^] neighbor (s) 
with equal strength. It is found that the CB depends on the range of coupling. In the 
global- coupling case, in which each ID map is coupled to all the other ID maps with equal 
strength, both the structure of the critical set and the CB are the same as those for the 
two-coupled case, irrespectively of A^. However, for the cases of nonglobal couplings of 
shorter range, a significant change in the structure of the critical set may or may not occur 
according as the coupling is linear or not. As examples of the linear and nonlinear nonglobal 
couplings, we study the linearly- and diffusively-coupled, nearest-neighbor coupling cases, 
respectively. For the linearly-coupled case, of the infinite number of period p-tupling routes 
for the global-coupling case, only the route ending at the zero-coupling critical point is left 
in the parameter plane. On the other hand, for the diffusively-coupled case, one critical line 
segment constitutes the critical set, as in the globally-coupled case. 

A. Stability of periodic orbits in many coupled maps 

Consider A^ symmetrically coupled ID maps with a periodic boundary condition, 

= F{Xmit),Xm+l{t), . . . 

m = l,...,A^, (4.1) 

36 



where is a positive integer larger than or equal to 2, x = (xi, . . . , xat) and cr is the cyclic 
permutation of x [i.e., crx = (x2, . . . ,Xi)]. Here Xm{t) is the state of the mth element at 
a discrete time t, and the periodic condition imposes Xm(t) = Xm+Nit) for all m. Like the 
two-coupled case with N = 2, the function F consists of two parts: 

F(x) = /(xi)+^?(x), (4.2) 

where / is an uncoupled ID map with a quadratic maximum at x = 0, and g is a coupling 
function. The unoupled ID map / satisfies the normalization condition ( |2.2| ), and the 
coupling function g obeys the condition 

g{x, . . . ,x) = for any x. (4.3) 

The A^-coupled map T has a cyclic permutation symmetry, 

CT-^Ta(x) = T(x) for all x, (4.4) 

where is the inverse of a. The set of all fixed points of a forms a symmetry line on 
which 

xi = ■ ■ ■ = xn- (4.5) 



It follow from Eq. ( [4.4| ) that the cyclic permutation a commutes with the map T, i.e., 
crT = Tcr. Hence the symmetry line becomes invariant under T, i.e., if a point x lies on the 
symmetry line, then its image T(x) also lies on it. An orbit is called a synchronous orbit if 
it lies on the symmetry line, i.e., it satisfies 

x^{t) = ■■■ = x^it) = x{t) for all t. (4.6) 

Otherwise, it is called an asynchronous orbit. Here we study only the synchronous orbits. 
They can be easily found from the uncoupled ID map, x{t + 1) = f{x{t)), because of the 
condition ([4.31). 

Consider an element, say the mth element, in the A^-coupled map T. Then the (m±5)th 
elements are called the 6th. neighbors of the mth element, where 1 < 5 < ^i^^^^) for even 
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(odd) A^. If the coupling extends to the i^th neighbor (s), then the number K is called the 
range of coupling. 

A general form of coupling for odd (A^ > 3) is given by 

K 



g{xi, ...,xn) 



2K + 



" m=—K 
K 



2K + 

K = h., 



- u{Xi+rn) - u{Xi) 



' m=—K 

N -1 



(4.7) 



where c is a coupling parameter and m is a function of one variable. Here the coupling 
extends to the Kth neighbors with equal coupling strength, and the function g satisfies the 
condition (|4.3| ). The extreme long-range interaction for K = ^^-^ is called a global coupling, 
for which the coupling function g becomes 



N 



m=l 



1 ^ 
m=l 



u{xi) 



(4.8) 



This is a kind of mean-field coupling, in which each element is coupled to all the other 
elements with equal coupling strength. All the other couplings with K < (e.g., nearest- 
neighbor coupling with K = 1) will be referred to as non-global couplings. The K = 1 case 
for A^ = 3 corresponds to both the global coupling and the nearest-neighbor coupling. 

We next consider the case of even A^ (A^ > 2). The form of coupling of Eq. ( [4.7|) holds 
for the cases of non-global couplings with K = 1, . . . , {N > 4). The global coupling for 



K = Y ^ 2) also has the form of Eq. 



but it cannot have the form of Eq. ( [4.7] ), 



because there exists only one farthest neighbor for i^" = unlike the case of odd A^. The 
K = 1 case for N = 2 also corresponds to the nearest-neighbor coupling as well as to the 
global coupling, like the A^ = 3 case. 

The stability analysis of an orbit in many coupled maps is conveniently carried out 
by Fourier-transforming with respect to the discrete space {m} pO|]. Consider an orbit 
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{xm{t) ', m = 1, . . . , A^} of the coupled maps ( |4.1|) . The discrete spatial Fourier transform 
of the orbit is: 

1 ^ 

m=l 

j = 0,l,...,N-l. (4.9) 

The Fourier transform ^^(t) satisfies ^j(t) = C,N~j{t) (* denotes complex conjugate), and the 
wavelength of a mode with index j is y for j < y and for j > 

To determine the stability of a synchronous orbit [xiit) = ■ ■ ■ = XN(t) = x{t) for all t], 
we consider an infinitesimal perturbation {6xm{t)} to the synchronous orbit, i.e., Xm{t) = 
x{t) + 6xm{t) for m = 1, . . . ,N. Linearizing the A^-coupled map (^4.1|) at the synchronous 
orbit, we obtain: 

N 

5x^{t + 1) = f\x{t))5x^{t) + Y.G^'\x{t)) 5x1+^-1, (4.10) 

1=1 

where 



^ d_g{xi,...,XN) 



(4.11) 

a;i=---=X]v=a; 



Hereafter the functions G'^'^'s will be called "reduced" coupling functions of g{xi, . . . ,xn)- 
Let S^j{t) be the Fourier transform of 6xmit), i.e., 

1 ^ 

m=l 

j=0,l,...,iV-l. (4.12) 

Here 6C,o is the synchronous- mode perturbation, and all the other S^j^s with nonzero in- 
dices j are the asynchronous-mode perturbations. Then the Fourier transform of Eq. (|4.10|) 
becomes: 

sut + I) = [f'ixit)) 

N 



1=1 



+ ^G'«(x(t))e2-('-i)^^/^]5e,(t), 

j = 0,l,...,iV-l. (4.13) 
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Note that all the modes S^/s become decoupled for the synchronous orbit. 

For a synchronous orbit with period q, its linear stability is determined by iterating the 
linearized map ( |4.13| ) q times: 



m=0 
N 

1=1 

j = 0,l,...,iV-l. (4.14) 
That is, the stability multipliers of the orbit are given by 

t=0 1=1 

j = 0,l,...,iV-l. (4.15) 

Here the first stability multiplier Aq is associated with the stability against the synchronous- 
mode perturbation, and hence it may be called the synchronous stability multiplier. On 
the other hand, all the other stability multipliers A/s (j 7^ 0) are called the asynchronous 
stability multipliers, because they are associated with the stability against the asynchronous- 
mode perturbations. 

A synchronous orbit becomes stable when it is stable against all the synchronous-mode 
and asynchronous-mode perturbations, i.e., the moduli of all stability multipliers are less 
than unity (|Aj| < 1 for j = 0, . . . , N — 1). Hence the stable region of the synchronous 
orbit in the parameter plane is bounded by the synchronous and asynchronous bifurcation 
lines determined by the equations Xj = ±1 for j = 0, . . . , — 1. When the Aq = 1(— 1) 
line is crossed, the synchronous orbit loses its stability via synchronous saddle-node (period- 
doubling) bifurcation. However, when the Aj (j 7^ 0) = 1(— 1) line is crossed, it becomes 
unstable via asynchronous pitchfork (period-doubling) bifurcation. 

It follows from Eq. (gj) that 

N 

J2G^'\x) = 0. (4.16) 
1=1 
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Hence the synchronous stabihty multipher Aq for j = becomes 

Ao = n/'(^W), (4.17) 
t=o 

which is just the stabihty multipher of the uncoupled ID map. While there is no coupling 
effect on Aq, the coupling generally affects asynchronous stability multipliers A/s of j ^ 0. 
The effect of the coupling on the asynchronous stability multipliers depends on the range of 
coupling, as wih be seen in the next two subsections. 



B. Global-coupling case 

In this subsection, we study the CB of period p-tuplings in many coupled maps with a 
global coupling. It is shown that both the structure of the critical set and the CB for the 
case of N globally-coupled maps are the same as those for the case of two coupled maps, 
independently of N . 

In case of the global coupling (^78|), the reduced coupling functions of Eq. ( |4.11| ) become: 



, [\-N)G{x) for/ = l, 
= { (4.18) 

G{x) for / ^ 1, 

where G{x) = jju'{x). Substituting 's into Eq. ( |4.15| ), we find that all the asynchronous 
stability multipliers are the same: 

9-1 

Ai = ■ ■ ■ = A^_i = n [fW)) - NG{x{t))] 
t=o 

= li[f'{4t))-cu'{xm. (4.19) 

i=0 

Hence there exist only two independent stability multipliers Aq and Ai (= Ai = • • • = Aat-i) 
for the global-coupling case. Note also that the values of Aq and Ai are independent of N and 
they are the same as those of two coupled maps. Thus the stability diagram of synchronous 
orbits of period {n = 0, 1, 2, . . .) in any globally coupled maps becomes the same as 
that of the two coupled maps. Consequently, the two parameter scaling factors associated 
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with scaling of the nonhnearity and couphng parameters are the same as those of the two 
coupled maps, independently of A^. That is, the CB of globally coupled maps becomes 
the same as that for the case of the two coupled maps, irrespectively of A^, which is also 
shown below by a renormalization analysis. 

We now follow the same procedure of Sec. [11 B| and straightforwardly extend the renor- 
malization results of the two coupled maps to many globally coupled maps. The rescaling 
operator B of Eq. ( p.lO|) becomes al for the case of N coupled maps, where I is the N x N 
identity matrix. Applying the period p-tupling operator Af of Eq. ( p^) to the A^-coupled 
maps (|4.1| ) n times, we obtain the n-times renormalized map T„ of the form. 



= fn{Xm{t)) +g it)), 

m = l,...,N. (4.20) 

Here /„ and Qn are the uncoupled and coupling parts of the n-times renormalized function 
Fn, respectively. They satisfy the following recurrence equations: 



/n+i(xi)=ail^)(^), 
a 

9n^i{^) = aF^^\-) - afi^\^ 
a a 



(4.21) 
(4.22) 



where 



flfK^l)= fn{f^'-'\x,)), 

Fjf\^)=FMP-'\^),F,!f-'\a^), . . . ,F(P-^)(a^-ix)), 



(4.23) 
(4.24) 



and the rescaling factor is chosen to preserve the normalization condition /n+i(0) = 1, i.e., 
= jif^TT^- Then, the recurrence relations ( f4.21| ) and ( [4. 221) define a renormalization 
operator TZ of transforming a pair of functions (/, g) : 



(4.25) 



(f ) 


= n 


(u 


\ 






^ 9n+l J 




^9n 


I 
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A critical map with the nonhnearity and couphng parameters set to their critical values 
is attracted to a fixed map T* under the iterations of the renormalization transformation 



T* : x^{t+ 1) = . . . ,x„_i(t)) 

= r{Xm{t)) + g*{Xra{t),Xm+l{t), • • • , 

m = l,...,iV, (4.26) 

where {f*,g*) is the fixed point of TZ with a = 1//*(1). Since /* is just the ID fixed map, 
only the equation for the coupling fixed function g* is left to be solved. 

As in the two-coupled case, we construct a tractable recurrence equation for the reduced 
coupling function G^''\x). That is, differenting the recurrence equation ( [4.25|) with respect 
to (Z = 2, . . . , A^) and setting xi = ■ ■ ■ = xn = x, we obtain 



(p),X. 



a 



a a a a 

N 



m=l 



1 = 2,. ..,N. 



(4.27) 



where F^fHx) = dF^^\^)/d: 



Xri 



[m 



1, . . . , A^). Note that these reduced cou- 



N 



pling functions satisfy the sum rule of Eq. (CT) (i-e., Y.G^}{x) = ^), and GW(x) 



G^^+^\x) [or equivalently, F^^)^ 



X] 



7(P) 
n,m+N 



1=1 



[x)] due to the periodic condition. 



In the global-coupling case, the initial reduced coupling functions {G^''\x)} satisfy 
Eq. ( [4.18| ), i.e., there exists only one independent reduced coupling function G{x). Then, 
it is easy to see that the successive images {Gll\x)} of {G^''\x)} under the transformation 



( ^ ) also satisfy Eq. ( ^ISj) , i.e. 



X) 



Gf)(x) = G„(x), 



[or equivalently, F^flix) 



7iP) 
n,N 



:x)^Fi^\x).] 



(4.28) 
(4.29) 



43 



Consequently, there remains only one recurrence equation for the independent reduced cou- 
pling function G{x): 

a 

a a a 

+f^'-'H-)Gn{fr'\-)). (4.30) 
a a 



Then, together with Eq. (|4.21| ), Eq. ( [4.30|) defines a reduced renormalization operator TZ of 



transforming a pair of functions (/, G) such that (/n+i, Gn+i) = T^ifn, Gn)- 

We look for fixed points of H such that = n{f*,G*). Here f*{x) is 

just the ID fixed function. Only the equation for G* is therefore left to be solved. Since 



the transformation ( |4.3CI| ) for G holds for any globally-coupled-map cases with > 2, it 
can be regarded as a generalized version of Eq. ( |2.21| ) for the two-coupled case. Comparing 
the expression in Eq. ( [4.30| ) with that in Eq. ( p.21| ), one can easily see that they are the 



same except for the factor A^. Making a change of the independent reduced coupling function 
G{x) jfG{x) [equivalents, F^p\x) ^F^p\x)], Eq. (|;3Up is transformed into Eq. {^^. 



Consequently, rescaling the solutions ( ^.331 ) for the two-coupled case with the scaling factor 
jj, one can obtain the solutions for the case of globally coupled maps: 

G*{x) = for all p, (4.31a) 

G*{x) = ^f*'{x) for all p, (4.31b) 

G*(x) = ^[r'(x)-1] for all p, (4.31c) 

G*{x) = ^ir'ix) + I] for odd p, (4.31d) 

G*{x) = ^f*'{x) for odd p. (4.31e) 

Thus, there exist three (five) fixed points (/*,(?*) of TZ for the case of even (odd) p, inde- 
pendently of N. 

For the same reason as for the two coupled maps [see Eq. ( p.36|) ], the critical stability 



multipliers have the values of the stability multipliers of the fixed point of the fixed map T*. 
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From Eqs. ( [4.17|) and ( [4.19| ), we obtain two independent critical stability multipliers Aq and 

A- 

K = r'{x), Xl = r'{x)~NG*{x), (4.32) 

where x is the fixed point of the ID fixed function [i.e., x = /*(£)] and Aq is just the critical 
stability multiplier A* of the uncoupled ID map. Substituting G*'s into Eq. ( [4.32| ), we obtain 
the same critical asynchronous stability multipliers A^'s as in the case of two coupled maps, 
as listed in Table ^. 

Consider an infinitesimal perutrbation to a fixed point {f*,G*) of the reduced 

renormalization operator TZ. Linearizing TZ at the fixed point, we obtain the recurrence 
equation for the evolution of {hn+i,^n+i) = '^ihn,^n)- As in the two-coupled case, 

the linearized operator C also has a semi-block form [i.e., hn+i{x) = [Cihn]{x) and <l'„+i(x) = 
[£2^n](a;) + [CshnKx)]. It follows from the reducibility of C into a semi-block form that one 
can find eigenvalues of Ci and C2 separately and then they give the whole spectrum of C 

All the fixed points (/*, G*) have a common relevant eigenvalue 6 of Ci (i.e., the relevant 
eigenvalue for the case of the uncoupled ID maps) associated with the critical scaling of 
the nonlinearity parameter of the uncoupled ID map. However, the relevant CE's of C2, 
associated with the critical scaling of the coupling parameter, depend on the kind of the 
fixed points, as in the case of the two coupled maps. Consider an infinitesimal perturbation 
$ to a fixed point G* of the recurrence equation ( [4.30|) . Linearizing Eq. ( [4.30|) at the fixed 



point G*, we obtain an equation for the evolution of $ [i.e., <l'„+i(x) = C2^n{x)]. Note again 
that Eq. ( [4.30| ) is transformed into the recurrence equation (|2.21|) for the two-coupled case 
under a mere scale change of the independent reduced coupling function G{x) — ^ j^G{x). 
Consequently, the CE equation [i.e., C2^*{x) = i^$*(x)] for the case of globally coupled 
maps becomes the same as that for the case of two coupled maps, independently of A^. Then, 
following the same procedure of Sec. |11 lj|, one can obtain the same CE's as those for the 



case of two coupled maps, as listed in Table |^. 
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C. Nonglobal-coupling cases 



In this subsection, we choose /(x) = 1 — Ax"^ as the uncoupled ID map and study the 
CB of the period p-tuphngs for the nonglobal-couphng cases. A significant change in the 
structure of the critical set may or may not occur according as the coupling is linear or not. 
As examples of the linear and nonlinear nonglobal couplings, we study the linearly- and 
diffusively-coupled, nearest-neighbor coupling cases, respectively. For the linearly-coupled 
case, only the zero-coupling critical point is left in the parameter plane, which is in contrast 
to the global-coupling case. On the other hand, for the diffusively-coupled case, one critical 
line segment constitutes the critical set, as in the global-coupling case. 

Consider a nonglobal coupling of the form ([4.7|) and define 



G(x) = — u'(x), (4.33) 



where 1 < K < ^^{^^y^) for even (odd) larger than 3. Then we have 



G'^'\x) 



-2KG{x) for / = 1, 
G{x) for 2 < / < 1 + or 

for A^ + 1 - a: < / < AT, 
otherwise. 



(4.34) 



Substituting G*^')'s into Eq. ( |4.15| ), we find that each stability multiplier \j (associated with 
the stabihty against the jth-mode perturbation) is given by 

A, = n[/' (^(^)) - S^{K,3)cu'{x{t))l (4.35) 

i=0 

where 

c (r. ■^- ■ 2^3k sin(2jr+l)f 

S]v(K,']) = > sm = 1 — ^ 4.36 

'-^^ 2K+lf^^ N {2K + l)smf ^ ' 

Hence, unlike the global-coupling case [see Eq. ( |4.19| )], the stability multipliers vary de- 



pending on the coupling range K as well as on the mode number j. Since SN{K,j] 
Siy{K, N — j), the stability multipliers satisfy 
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A, = A^_,-, J =0,l,...,iV-l. (4.37) 

Thus it is sufficient to consider only the case of < j < y {^^^) even (odd) N. 
Comparing the expression in Eq. ( [4. 35] ) with that in Eq. ( ^.191 ) for j ^ 0, one can easily 



see that they are the same except for the factor Si\f{K,j). Consequently, making a change 
of the coupling parameter c Sm{k j) ' stability multiplier Xj for the nonglobal-coupling 
case of range K becomes the same as that for the global-coupling case. 

For each mode with nonzero index j, we consider a region in the parameter plane, in 
which a synchronous orbit is stable against the perturbations of both modes with indices 
and j. This stable region is bounded by four bifurcation curves determined by the equations 
Ao = ±1 and Xj = ±1, and it will be denoted by Un- For the case of global coupling, those 
stable regions coincide, irrespectively of and j, because all the asynchronous stability 
multipliers A^'s (j ^ 0) are the same, independently of N. The stable region for this global- 
coupling case will be denoted by Ug- Note that Uq itself is just the stability region of the 
synchronous orbit, irrespectively of A^, because the synchronous orbit is stable against the 
perturbations of all the synchronous and asynchronous modes in the region Ug- 

However, the stable regions Un^s vary depending on the coupling range K and the mode 
number j for the nonglobal-coupling cases, i.e., Un = UN{K,j). To find the stability region 
of a synchronous orbit in an A^-coupled map with a given K, one may start with the stability 
region Ug for the global-coupling case. Rescaling the couphng parameter c by a scaling factor 
s^l^x j) Gach nonzero j, the stable region Ug is transformed into a stable region UN{K,j). 
Then the stability region of the synchronous orbit is given by the intersection of all such 
stable regions U^^s. A significant change in the stability diagram of the synchronous orbits 
of period {n = 0, 1, 2, . . .) may or may not occur according as the coupling is linear or 
not, as will be seen below. 

As the first example, we study the linearly-coupled, nearest-neighbor coupling case with 
K = 1, in which the coupling function is 

g{xi, . . . , Xn) = -ix2 + Xn - 2xi) for A^ > 3. (4.38) 

o 
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This kind of coupling can be regarded as a generalized version of Eq. ( p.l|) for the linearly- 
coupled case for N = 2. (Note that as mentioned above, the cases of = 2 and 3 correspond 
to the global-coupling case.) For K = 1, the scaling factor Sj^{k j) '^^ ^'i- ( |4-36| ) becomes 

SN{l,j) = lsm\^). (4.39) 

This scaling factor ^^^^ has its minimum value, \{ ^^^^2 jl_ ) at jmin = ^{^^^) for even 
(odd) A^. We also note that the minimum value for odd N depends on N , but d& N oo 
it converges to the minimum value for the case of even A^. 

Rescaling the coupling parameter c with the scaling factor -s^^^TJ)i the stability region 
Ug for the N = 2 and 3 cases of global coupling is transformed into the region UnO-iJ) for 

> 3. Then, as a result of the intersection of all such regions ?77v(l, j)'s, only the region 
t/iv(l, jmin) including a c = line segment is left as the stability region of a synchronous 
orbit. Consequently, of the infinite number of period p-tupling routes for the global-coupling 
case, only the Zp route ending at the zero-coupling critical point (A^^^ 0) remains. Thus 
only the zero-coupling point is left as a critical point in the parameter plane. An example 



for the period-tripling case with p = 3 is shown in Fig. The largest stability region 
corresponds to the stability region for the N = 2 and 3 cases of the global coupling, while 
the smallest one corresponds to that for the case of even A^ (A^ > 2). Between them, there 
exist stability regions for the case of odd A^ (A^ > 3) (e.g., see the A^ = 5 case in Fig. ^1]) . 

We now examine the CB near the zero-coupling critical point for the case of the linearly- 
coupled, nearest-neighbor coupling. Consider a self-similar sequence of parameters (yl„,c„), 
at which the synchronous orbit of level n has some given stability multipliers, in the Zp route 
for the global-coupling case. Rescaling the coupling parameter with the factor „ — r, this 
sequence is transformed into a self-similar one for the case of the linearly-coupled, nearest- 
neighbor coupling. Thus the "width" of each stability region in the Zp route for the case of 
the global coupling is reduced to that for the case of the linearly-coupled, nearest-neighbor 
coupling by the scaling factor g^^^ij y while the "heights" of all stable regions in the Zp 



route remain unchanged [0. It is therefore obvious that the critical scaling behavior near 
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the zero-coupling critical point for the case of the linearly-coupled, nearest-neighbor coupling 
is the same as that for the global coupling case. That is, the height and width hn and Wn of 
the stability region of level n geometrically contract in the limit of large n, 

hn ~ Wn ~ for large n, (4.40) 

where S and a are the scaling factors of the nonlinearity and coupling parameters, respec- 
tively. As an example, see again Fig. ^ and note that Figs, ^(a), ^(b) and |2l|(c) nearly 



coincide near the zero-coupling critical point except for small numerical differences. 

The results of the linearly-coupled, nearest-neighbor coupling with K = 1 extend to all 
the other linearly-coupled, nonglobal-coupling cases with 1 < K < ^i^^T^) even (odd) 
A^. For each nonglobal-coupling case with K > 1, we first consider a mode with index j^^^ 
for which the scaling factor g^^j^ j-^ becomes the smallest one. Here the value of jmin varies 
depending on the range K. Like the K = 1 case, only the region f/Af(A', jmin) including a 
c = line segemnt is left as the stability region of a synchronous orbit. Consequently, only 
the zero-coupling point remains as a critical point in the parameter plane, and the CB near 
the zero-coupling critical point is also the same as that for the global-coupling case. 

Finally, as an example of the nonlinear coupling, we study the diffusively- coupled, 
nearest-neighbor coupling case with = 1, in which the coupling function is given by 

gix,, ...,x^) = ^[f{x2) + fM - 2/(xi)] for N > 3. (4.41) 

This kind of coupling can be regarded as a generalized version of Eq. ( |3.32| ) for the 
dissipatively-coupled case for N = 2. For this diffusively-coupled case, the stability multi- 
pliers of Eq. (|4.35| ) for a synchronous g-periodic orbit of level n {q = p") become 

A, = [l-c5^(l,j)]^Ao. (4.42) 

Like the linearly-coupled case, rescaling the coupling parameter c with the minimum 
scaling factor . y each stability region Ug of level n for the global-coupling case is 
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transformed into the stability region [i.e., the region Un{1, jmin)] of level n for the diffusively- 
coupled case. The stability region U]\f{l,jmm) is bounded by four bifurcation curves deter- 
mined by the equations = ±1 for j = 0, jmin- An infinite sequence of such stability regions 
converges to a critical line segment joining two points cf (= 0)and c* (= ) ) on the 

A = A'^p} line. Consequently, one critical line segment constitutes the critical set, like the 
global-coupling case. An example for the period-tripling case for even (A^ > 2) is shown 
in Fig. 

As shown above, the stability diagram for the diffusively- coup led case is essentially the 
same as that for the globally-coupled case, except for the scale in the coupling parameter 
c. Hence the critical scaling behavior of the nonlinearity and coupling parameters and so 
on becomes the same as that for the global-coupling case (for details on the CB, refer to 
Sec. [Ill Bp . Finally, we briefly discuss the critical asynchronous stability multipliers. At the 
zero-coupling critical point (A^^O) and interior critical points, they are the same as those 
for the global-coupling case (i.e., = ■ ■ ■ = A^.^ = Aq at {A^^, 0), and A^ = ■ ■ ■ = A^_i = 
at interior critical points.). However, at the right end (A^^c*), the critical asynchronous 
stability multiplier for j = jmin is given by A*^^^.^^ = Aq (even q) and — Aq (odd q), but all other 
A*'s (j 7^ jmin) are zero, which is somewhat different from that for the global-coupling case 
[i.e., XI = ■ ■ ■ = = Aq (even q) and — Aq (odd q)]. Like the linearly-coupled case, the 

results for the diffusively-coupled, nearest-neighbor coupling case with K = 1 can be also 
extended to all the other diffusively-coupled, nonglobal-coupling cases with 1 < K < ^i^^) 
for even (odd) A^. 

V. SUMMARY 

The CB of all period p-tuplings {p = 2,3,4,...) are studied in (A^ = 2,3,4,...) 
symmetrically coupled ID maps. The two-coupled case with A^ = 2 is first investigated by 
a renormalization method. We find three (five) kinds of fixed points of the renormalization 
operator and their relevant CE's associated with coupling perturbations for the case of 
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even (odd) p. We next consider two kinds of couplings, linear and nonlinear couplings. As 
examples of the linear- and nonlinear-coupling cases, we study the linearly- and dissipatively- 
coupled maps, respectively, and confirm the renormalization results. The structure of the 
critical set varies depending on the nature of coupling. For the linearly-coupled case, an 
infinite number of the critical line segments, together with the zero-coupling critical point, 
constitute the critical set, while for the dissipatively-coupled case, the critical set consists of 
the only one critical line segemnt. The CB also depends on the position on the critical set. 
For even (odd) p, three (four) kinds of fixed points govern the CB for the linearly-coupled 
case, while only two (three) kinds of fixed points govern the CB for the dissipatively-coupled 
case. Finally, the results of the two coupled maps are extended to many coupled maps with 
N > 3, in which the CB depends on the range of coupling. In the global-couphng case, 
both the structure of the critical set and the CB are the same as those of the two-coupled 
case, independently of N. However, for the nonglobal-coupling case, a significant change in 
the structure of the critical set may or may not occur according as the coupling is linear 
or not. For the linearly-coupled case, of the infinite number of period p-tupling routes for 
the global-coupling case, only the route ending at the zero-coupling critical point is left in 
the parameter plane. On the other hand, for the diffusively-coupled case, one critical line 
segment constitutes the critical set, as in the globally-coupled case. 
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TABLES 

TABLE I. Reduced coupling fixed functions G*{x), relevant CE's u and the critical asyn- 
chronous stability multipliers A| in all the period p-tupling cases are shown for the case of two 
coupled maps. In the second column, "e" ("o") denotes the even (odd) period-p tuplings. The first 
three ones for G* {x) exist for all p, while the last two ones exist only for odd p. Note also that 

the case G*{x) = ^f*'{x) has no relevant CE's, and a and A* are the orbital scaling factor and the 



critical stability multiplier for the ID case, respectively. 

G*{x) p u XI 

e, o a, p X* 

I /*' e, o nonexistent 

Ur'{x)-1] e,o p 1 

i[r'(x) + i] o p -1 

r'(x) o a.p -X* 



TABLE IL In the Z3 route, we followed a sequence of parameters {An, Cn) at which the pair 
of stability multipliers {Xo,n,^i,n) of the orbit of period 3" is (—1, 1). This sequence converges to 
the zero-coupling critical point with the scaling factors shown in the second and third 

columns. 



n 6n fJ-n 

3 55.264 789 71 -9.272 61 

4 55.245 77151 -9.279 20 

5 55.247110 93 -9.276 71 

6 55.247 020 84 -9.277 55 

7 55.247 026 98 -9.277 27 

8 55.247 026 56 -9.277 36 

9 55.247 026 59 -9.277 33 
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TABLE III. For the case of the Zj, route, scahng factors and /X2,n in the two-term scaUng 
for the couphng parameter are shown in the second and third columns, respectively. A product of 
them, — is shown in the fourth column. 



4 -9.277 396 06 24.578 89 3.50179 

5 -9.27733731 27.768 52 3.099 52 

6 -9.27734136 28.502 38 3.019 72 

7 -9.277 341 10 28.650 77 3.004 08 

8 -9.277 34112 28.68153 3.000 85 



TABLE IV. In the period-tripling case, scaling factors ri^„ and r2,n in the two-term scaling for 
the slope Sn of the asynchronous stability multiplier Ai^„ at the zero-coupling critical point are 
shown in the second and third columns, respectively. 



■n ri,„ r2,n 

4 -9.277335 5434 2.927 8 

5 -9.2773415010 2.984 5 

6 -9.277341089 3 2.996 7 

7 -9.2773411174 2.999 3 

8 -9.277341116 8 2.999 5 
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TABLE V. We followed, in the leftmost L3 route, two self-similar sequences of parameters. 
One sequence of parameters {An,Cn), at which Ao,n = — 1 and Ai,„ = 0.8, converges to the left 
end {Aoo ,ci) of the critical line segment with the scaling factor fin of the coupling parameter 
shown in the second column. The other sequence of parameters (A„,c„), at which Ao,n = — 1 
and Xi^n = —0.8, converges to the right end {A^ ,Cr) with the scaling factor /x„ of the coupling 
parameter shown in the fourth column. In both cases the scaling factors are the same. The 
convergence of the sequence {Ai^^} to its limit values A* at the left (right) end is also shown in the 
third (liflli) cohuuu. Xotc; llial the values of AJ"s at IkjiIi cuds are dillcrcnl. 



3 3.664 5 0.996 763 57 1.836 2 -1.010 285 72 

4 2.946 7 1.000 44409 3.0573 -0.998 599423 

5 2.990 8 0.999 93919 2.968 2 -1.00019198 

6 2.9941 1.000 000 33 2.9970 -0.999 973 71 

7 2.998 3 0.999 998 86 2.998 -1.000 003 60 

8 2.999 4 1.000 00016 2.999 5 -0.999 999 51 

9 2.999 8 0.999 999 98 2.999 8 -1.000 000 07 
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TABLE VI. In the Z4 route, we followed a sequence of parameters {An, Cn) at which the pair 
of stability multipliers (Ao,n! Ai,n) of the orbit of period 4" is (—1, 1). This sequence converges to 
the zero-coupling critical point {J^^ , 0) with the scaling factors shown in the second and third 
columns. 



n (5„ \in 

2 981.297 281 -38.84191 

3 981.596 790 -38.81617 

4 981.594 965 -38.819 38 

5 981.594 977 -38.819 04 

6 981.594 976 -38.819 08 



TABLE VII. For the case of the ^4 route, scaling factors /ii.n and //2,n in the two-term scaling 
for the coupling parameter are shown in the second and third columns, respectively. A product of 
them, is shown in the fourth column. 

n llx,n [ll,n ^ 



3 -38.819 02199 311.92 4.83 

4 -38.819 074 56 372.98 4.04 

5 -38.819 074 24 377.22 3.99 



TABLE VIII. In the period-quadruphng case, scaling factors t\ji and r2 n in the two-term 
scaling for the slope Sn of the asynchronous stability multiplier Ai^„ at the zero-coupling critical 

point arc shown in the second and third columns. rcs]")(xMi\'cly. 



n ri^n r2,n 

2 -38.819 023 65 3.593 

3 -38.819 07465 3.972 

4 -38.819 074 29 3.998 
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TABLE IX. We followed, in the leftmost L4 route, two self-similar sequences of parameters 
{An,Cn), at which Ao,n = —1 and Ai,„ = 0.8. One sequence converges to the left end {A^,ci) of 
the critical line segment with the scaling factor of the coupling parameter shown in the second 
column. The other sequence converges to the right end {A^,Cr) with the scaling factor /x„ of the 
coupling parameter shown in the fourth column. In both cases the scaling factors are the same. 
The convergence of the sequence {Ai^n} to its limit values AJ at the left (right) end is also shown 
in the third (fifth) column. Note that the values of A^'s at both ends are the same. 



n l^n Ai,„ IJ,n Xl,n 

2 5.013 1.010 032 65 50.70 0.987196 90 

3 4.051 0.999 630 89 3.674 1.000476 66 

4 3.963 1.000 013 65 3.978 0.999 982 38 

5 3.992 0.999 999 50 3.993 1.000 000 65 

6 3.998 1.000 000 02 3.998 0.999 999 98 



TABLE X. Independent reduced coupling fixed functions G*{x), relevant CE's v and indepen- 
dent critical asynchronous stability multipliers A* in all the period p-tupling cases are shown for 
the case of N globally coupled maps. The first three ones for G*{x) exist for all p, whereas the last 
two ones exist only for odd p. Here a and A* are the orbital scaling factor and the critical stability 
multiplier for the ID case, respectively. Note that u and A^ for each G*{x) are also the same as 



those for the two-coupled case, independently of N. 

G*{x) V \\ 

a, p A* 

f*' nonexistent 

^f*'{x) a, p -A* 
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FIGURES 

FIG. 1. Stability diagram of the synchronous 3"-periodic orbits of the lowest two levels n = 0,1 
in two linearly coupled maps. Each periodic orbit of level n is born via its own saddle-node bifur- 
cation. Its stable regions, denoted by g = 3", are bounded by four bifurcation curves determined 
by Xi = ±1 for i = 0,1. The horizontal and non-horizontal solid (short-dashed) boundary lines 
(i.e., the Aq = —1 (1) and Ai = —1 (1) bifurcation curves) correspond to the synchronous and 
asynchronous period-doubling (synchronous saddle-node and asynchronous pitchfork) bifurcation 
curves, respectively. 

FIG. 2. Stability diagram of the synchronous 3"-periodic (n = 1, 2, 3) orbits of level n in two 
linearly coupled maps. Each periodic orbit of level n is born via its own saddle-node bifurcation. Its 
stable regions, denoted by g = 3", are bounded by four bifurcation curves determined by Aj = ±1 
for i = 0,1. The solid and short-dashed boundary lines represent the same as those in Fig. 1. The 
stability diagrams starting from the central, right, and left stability regions of level 1 are shown in 
(a), (b), and (c), respectively. See the text for other details. 

FIG. 3. Stability regions of the synchronous orbits of period g = 3" (n = 1, 2, 3) in two linearly 
coupled maps. The cases n = 1,2, and 3 are shown in (a), (b), and (c), respectively. Its stable 
regions, denoted by g = 3", are bounded by four bifurcation curves determined by Aj = ±1 for 
i = 0, 1. The solid and short-dashed boundary lines represent the same as those in Fig. 1. The 
scaling factors used in (b) and (c) are 5 = 55.247026 and a = -9.277341. 

FIG. 4. Plots of the asynchronous stability multipliers Ai^„(^to,c) versus c near the 
zero-coupling critical point for the synchronous orbits of period g = 3" (n = 2, 3, 4). 

FIG. 5. Plots of the asynchronous stability multipliers Xi,n{A'^,c) versus c near the leftmost 
critical line for the synchronous orbits of period g = 3" (n = 2, 3, 4). 
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FIG. 6. Plot of the asynchronous Lyapunov exponent ai{A^',c) versus c near the leftmost 
critical line. The plot consists of 200 c values, each of which is obtained by iterating the map 
100 000 times to eliminate transients and then averaging over another 500 000 iterations. 

FIG. 7. Attractors near the leftmost critical line: (a) a synchronous attractor inside the critical 
line and asynchronous attractors outside the critical line for (b) c = —3.590 5 and (c) c = —3.4824. 
For each case, the map is iterated 100 000 times to eliminate transients and the next 10 000 iterations 
are plotted. 

(3) 

FIG. 8. Plots of the asynchronous Lyapunov exponent ai{Aoo ,c) versus c near the (a) left 
and (b) right ends of the leftmost critical line. Each plot consists of 50 £ values, each of which 
is obtained by iterating the map 100 000 times to eliminate transients and then averaging over 
another 500000 iterations. 

FIG. 9. Stability diagram of the synchronous 4"-periodic orbits of the lowest two levels n = 0, 1 
in two linearly coupled maps. Each periodic orbit of level n is born via its own saddle-node bifur- 
cation. Its stable regions, denoted by ^ = 4", are bounded by four bifurcation curves determined 
by Aj = ±1 for i = 0, 1. The solid and short-dashed boundary lines represent the same as those in 
Fig. 1. 

FIG. 10. Stability diagram of the synchronous 4"-periodic orbits of level n (n = 1, 2, 3) in two 
linearly coupled maps. Each periodic orbit of level n is born via its own saddle-node bifurcation. Its 
stable regions, denoted hy q = 4", are bounded by four bifurcation curves determined by Aj = ±1 
for i = 0,1. The solid and short-dashed boundary lines represent the same as those in Fig. 1. The 
stability diagrams starting from the central, right, and left stability regions of level 1 are shown in 
(a), (b), and (c), respectively. See the text for other details. 
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FIG. 11. Stability regions of the synchronous orbits of period g = 4" (n = 1, 2, 3) in two hnearly 
coupled maps. The cases n = 1,2, and 3 are shown in (a), (b), and (c), respectively. Its stable 
regions, denoted by g = 4", are bounded by four bifurcation curves determined by Aj = ±1 for 
i = 0, 1. The solid and short-dashed boundary lines represent the same as those in Fig. 1. The 
scaling factors used in (b) and (c) are 5 = 981.594976 and a = -38.819 074. 

FIG. 12. Plots of the asynchronous stability multipliers \i,n{A^£ ,c) versus c near the 
zero-coupling critical point for the synchronous orbits of period 4" (n = 1, 2, 3). 

FIG. 13. Plots of the asynchronous stability multipliers Ai,„(A^\c) versus c near the leftmost 
critical line for the synchronous orbits of period 4" (n = 2, 3, 4). 

FIG. 14. Plot of the asynchronous Lyapunov exponent ai{A^(x!,c) versus c near the leftmost 
critical line. The plot consists of 200 c values, each of which is obtained by iterating the map 
100 000 times to eliminate transients and then averaging over another 500 000 iterations. 

FIG. 15. Attractors near the leftmost critical line: (a) a synchronous attractor inside the 
critical line and asynchronous attractors outside the critical line for (b) c = —3.877059 and (c) 
c = —3.888 059. For each case, the map is iterated 100 000 times to eliminate transients and the 
next 10 000 iterations are plotted. 

FIG. 16. Plots of the asynchronous Lyapunov exponent ai{A^,c) versus c near the (a) left 
and (b) right ends of the leftmost critical line. Each plot consists of 50 e vahies, each of which 
is obtained by iterating the map 100000 times to eliminate transients and then averaging over 
another 500000 iterations. 
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FIG. 17. Stability diagram of the synchronous orbits for the (a) period-triphng (p = 3) and (b) 
period-quadruphng (p = 4) cases in two dissipatively coupled maps. Each periodic orbit of level n 
is born via its own saddle-node bifurcation. Its stable regions, denoted hy q = (n = 1, 2, 3), are 
bounded by four bifurcation curves determined by Aj = ±1 for z = 0, 1. The solid and short-dashed 
boundary lines represent the same as those in Fig. 1. 

FIG. 18. Fixed points of the recurrence equations (a) (3.35) and (b) (3.36). The intersection 
points between the curves of Eqs. (3.35) and (3.36) and the line Cn = Cn-i are just the fixed points, 
denoted by solid circles. 

FIG. 19. Plots of the asynchronous stability multipliers Ai,„(A^\c) versus c for the syn- 
chronous orbits of period p'^ (n = 1,2,3) near the critical line for the (a) period-tripling {p = 3) 
and (b) period-quadrupling (p = 4) cases in two dissipatively coupled maps. 

FIG. 20. Asynchronous attractors outside the critical line for (a) c = —0.001 5 in the pe- 
riod-tripling case and (b) c = —0.000 1 for the period-quadrupling case. For each case, the map is 
iterated 100 000 times to eliminate transients and the next 10 000 iterations are plotted. 

FIG. 21. Stable regions of the synchronous 3'*-periodic orbits of level n = 1,2, 3 in iV 
{N = 2,3,4,...) coupled ID maps for the linearly-coupled, nearest-neighbor coupling case with 
K = 1. The cases of n = 1, 2, and 3 are shown in (a), (b), and (c), respectively. The symbols ^J"^, 
Qj^ , and Qj'^ with j = and jmin represent the period-doubling, saddle-node, and pitchfork bifur- 
cation lines associated with the stability of a synchronous orbit of period q against the perturbation 
of the jth mode. The scaling factors used in (b) and (c) are S = 55.247026 and a = —9.277341. 
For other details, see the text. 
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FIG. 22. Stable regions of the synchronous 3"-periodic orbits of level n = 1, 2, 3 in even- AT 
{N > 2) coupled ID maps for the diffusively-coupled, nearest- neighbor coupling case with K = 1. 
The symbols qj'^ , Qj'^, and gj"^ with j = and jjv represent the period-doubling, saddle-node, 
and pitchfork bifurcation lines associated with the stability of a synchronous orbit of period q 
against the perturbation of the jth mode. 
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